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INTEGRAL CONVEX POLYHEDRA 
AND 


AN APPROACH TO INTEGRALIZATION* 


Abstract 


Many combinatorial optimization problems may be formulated as 
integer linear programming problems -- that is, problems of the form: 
given a convex polyhedron P contained in the non-negative orthant 
of n-dimensional space, find an integer point in P which maximizes 
(or minimizes) a given linear objective function. Well known 
linear programming methods would suffice to solve such a problem if: 


(i) P is an integral convex polyhedron, or 


(ii) P is transformed into the integral convex polyhedron 
that is the convex hull of the set of integer points 
in P, a process which is called integralization. 


This thesis provides some theoretical results concerning integral 
convex polyhedra and the process of integralization. Necessary and 
sufficient conditions for a convex polyhedron P to have the integral 
property are derived in terms of the system of linear inequalities 
defining P. A number-theoretic method for integralizing two-dimensional 
convex polyhedra is developed which makes use of a generalization of 
the division theorem for integers. The method is applicable to a 
restricted class of higher dimensional polyhedra as well. 


*This report reproduces a thesis of the same title submitted to the 
Electrical Engineering Department, Massachusetts Institute of Technology, 
in partial fulfillment of the requirements for the degree of 

Doctor of Philosophy, August 1970. 
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BIOGRAPHICAL NOTE 


CHAPTER 1 


INTRODUCTION 


It is well known that many optimization problems of a combinatorial 
nature may be formulated as integer programming problems. Several authors 
have written on this subject, including Balinski [1] and [2], Dantzig [5], 
Hoffman [14], and Hu [17]. Outstanding examples are the problems of 
constructing a maximal flow in a transport network and constructing a 
maximal matching in a bipartite graph. In fact, these two examples lead 
very naturally to integer linear programming formulations which may be 
solved by standard linear programming methods. Other combinatorial 
optimization problems, such as the traveling salesman problem or the 
problem of properly coloring the regions of a map with a minimum number 
of different colors, lead to integer linear programming formulations 
which cannot be solved by standard linear programming methods. They 
require different and more complex solution methods. 

In 1958, Gomory [8] developed the first general algorithm for 
solving integer linear programming problems. There has been considerable 
interest in integer linear programming since Gomory's pioneering work, 
and a variety of methods for dealing with the problem have been 
developed. Balinski [1] and [2] presents an extensive survey of many 


of these methods and describes some computational experience with them. 


In 1965, Gomory [9] brought a new algebraic approach to bear on 
the problem. He showed that when certain simplifying assumptions are 
made, the integer linear programming problem can be transformed into 
an optimization problem involving a finite Abelian group. The 
algebraic ideas contained in his approach have been further developed 
and applied by Gomory [10] and [11], Shapiro [23] and [24], Hu [16], 
and others. 

Although a good deal of success in dealing with integer linear 
programming problems has been achieved, we believe there remains 
much to be learned in this area. Motivated primarily by the 
importance of integer linear programming as a tool for solving 
difficult combinatorial optimization problems, this thesis seeks 
to provide some relevant theoretical results which may contribute 
to a better understanding of the problem. 

The thesis is organized into seven chapters. In the remainder 
of this first chapter we review some n-dimensional geometric concepts 
and definitions related to convex polyhedra. We then give some 
graphical examples to illustrate the concept of an integral convex 
polyhedron and the concept of transforming a convex polyhedron into 
an integral convex polyhedron — a process which we call integralization. 
We conclude the chapter with a discussion of the relevance of integral 
convex polyhedra and the process of integralization to the integer 
linear programming problem, and a description of some of the related 


research of others. 


In Chapter 2 we derive several results concerning convex polyhedra 
having the integral property. In particular, we present a necessary 
and sufficient condition for a system of linear equations with integer 
coefficients and integer right hand sides to have an integer solution. 
We then use this to extend the results of Hoffman and Kruskal [15] 
and obtain a necessary and sufficient condition for a particular 
convex polyhedron to have the integral property. 

In Chapter 3 we define a new kind of integral property which we 
call the totally integral property. This definition is motivated by 
two considerations. First, the class of totally integral polyhedra 
is a proper subclass of the class of integral polyhedra. Second, 
the class of linear programming polyhedra having the totally integral 
property can be characterized algebraically, by applying some of the 
group theoretic ideas developed by Gomory [9]. 

In the next three chapters we turn to the process of integralization. 
We shall classify polyhedra using two parameters n and k. By an 
(n,k)-polyhedron we shall mean an n-dimensional polyhedron defined 
by a system of linear inequalities whose coefficient matrix has rank k, 
where k < n. 

In Chapter 4 we show that an (n,k)-polyhedron, where n is 
arbitrary, n =k, can always be transformed into a (k,k)-polyhedron 
by a transformation which preserves integralization. This allows us to 
focus on the problem of integralizing (k,k)-polyhedra without loss 


of generality. As we shall see, if we know how to integralize a 


(k,k)-polyhedron then we know how to integralize an (n,k)-polyhedron 
for arbitrary n >k. We then discuss the simple problem of 
transforming and integralizing (n,1)-polyhedra. We show that this 
is accomplished by using the greatest integer function applied to 
rational numbers, an application of the division theorem for integers. 
In Chapters 5 and 6 we deal with integralization of (2,2)-polyhedra. 
We introduce a generalization of the division theorem for integers 
which applies to ordered pairs of integers and use this, together 
with an associated division process, to integralize such polyhedra. 
Chapter 7 concludes the thesis with a discussion of the difficulty 
we encountered in attempting to generalize the integralization process 
to make it apply to (k,k)-polyhedra, for k => 3. We then offer some 


suggestions for future research in this area. 


1.1 Geometric Concepts and Definitions 


In this section we review some n-dimensional geometric concepts 
and definitions. Most of the material is well known, with the exception 
of the subsection on integral convex sets, and may be found in one or 
more of the following references: Benson [3], Bonnice and Klee [4], 
Goldman [6], Goldman and Tucker [7], Grunbaum [12], Hadley [13], 


Hoffman and Kruskal [15], Klee [19], and Weyl [26]. 


Let R be the real field. Let R” be the n-dimensional vector space 
over R. We refer to x € ae x = (x) > Xoo sees x)» as a real point 


or real n-vector (column). 


1.1.1 Dimension, Boundedness, and Convexity 


Oo il k ‘ 
Let x , X , «+e, X be real n-vectors. A real n-vector x is 
‘ ‘ ‘ ; Oo ol k 
said to be a linear combination of x , x, ..., x if 


0 1 k 
x = Xo x + Ny K tb eee + Me x 


for some real numbers Nop Ny? pies 5S If the A's satisfy the 


Kk 


constraint, 


Ag tra, toe tA al 


, : . ‘ ‘ Oo 1 k 
then x is said to be an affine combination of x , K , «.«, xX» 


If the X's satisfy the additional constraint, 


Ne 20, for i = 0, 1,...,k 


, j . j 0 1 k 
then x is said to be a convex combination of kK , XK , «ee, X- 


A set X = eae es waa -) is said to be linearly dependent 
if 0, the all zero n-vector, is a linear combination of oe sere eeey Pi 


in which some ve #0. If, in addition, 


+eee tr = O 
Xo + My k 


10 


then X is said to be affinely dependent. If X fails to be linearly 
(affinely) dependent then X is said to be linearly (affinely) 
independent. We shall say that X is a dependent (independent) 

set iff X is affinely dependent (independent). It can be 

shown that X is a dependent (independent) set iff the set 


0 2 0 0 
Geog 5 MERA. wey es } is linearly dependent (independent). 


3 : 
Example 1-1 Independent sets in R are: a single real point, 
two distinct real points, three real points not contained in a 
single line, four real points not contained in a single plane. oO 
n 
For the following definitions we let X be any subset of R , 
either a finite subset or an infinite subset. 
X is k-dimensional, where 0 < k < n, if X contains an independent 
set of k+l real points but none of k+2 real points. 
é : : 7 + 
X is bounded if there exist two real n-vectors x and x_ such 
- + : F , 
that, for all x € X, x < * < x , where < is componentwise. Otherwise 
X is unbounded. 
: F 1 2 
X is convex if, for all x", x ¢€ X and for all real dr, Oe A<l, 
1 2 n. : 
uN x + (1-\) x € X. In words, a convex subset X ¢ R is one which 
contains the line segment joining any two of its points. 
A s n. A 
An intersection of convex subsets of R is again a convex subset 


of R. 
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2 
Example 1-2 In R', any single real point is a O-dimensional bounded 
2 
convex set. R itself is a 2-dimensional unbounded convex set. Fig. i-l 
shows several subsets of R- which have various combinations of these 


properties. 


Bounded Unbounded 


1-dimensional 


2-dimensional 


Fig. 1-1 
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G ibe Bee Convex Hull 


Given an arbitrary subset X ¢ Re which may or may not be convex, 
it is always possible to find a convex subset of R which contains X. 
R” itself will do the job. The convex hull or convex closure of X, 
denoted by H(X), is defined to be the intersection of all convex 
subsets of R" which contains X. H(X) is unique. Intuitively, H(X) 
is the smallest convex subset of R” which contains X. If X is 


itself convex, then H(X) =X. 


2 
Example 1-3 Fig. 1-2 illustrates several subsets of R= and, 


immediately to the right of each, their convex hulls. 5 


(a) 7 
oe, 
7 (d) 
(b) 
. ; (e) 
(ec) 


Fig. 1-2 
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If X is a finite subset of R", X = ae i wee, x 3, then it can 
be shown (see Benson [3]) that the convex hull of X is given by, 
m m 
= n 2 i iss a 
HQ) S| mek fee) eK, Oye, ya, +1} (1-1) 
i= i=1 


where Nos Ay» eetig mn are real numbers. In words, (1-1) says that 
H(X) is the set of all convex combinations of X. 

Another characterization of the convex hull of a subset of R", 
similar to (1-1) but more general in that it applies to infinite 
subsets as well, is given by Caratheodory's Theorem (see Bonnice 


and Klee [4]). Given a subset X ¢ R", the theorem states that 


n 
AG) {een fe =: ). Rea hy SO Rey Seay GRD) 


bs 1, 
i=0 i=0 


where Ao» Mae Seeks nO are real numbers. In words, (1-2) states that 
H(X) is the set of all convex combinations of all finite subsets of 


X containing n+l or fewer points. 


1.1.3 Flats 


1 


0 
Lets Ke (ey ey aes x} be an independent set of real points 


in p. where 0 < k < n. A k-flat 7 in R is a k-dimensional convex 


subset of R", 
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k k 


m= {xen yx= ) a,x, Ya, =1} i535 
i=0 i=0 


where Nos Mie oe Me are real numbers. + is the set of all affine 
combinations of X. X is said to generate 7. A k-flat in R” is also 


known as a k-dimensional affine subspace of Ro 


Equivalently, a k-flat 7 is a set, 
n 
n={xeR | A x =b} (1-4) 


where A is a fixed real (n-k) y n matrix of rank n-k, and b is a 
fixed real (n-k)-vector. Both x and b are column vectors, and A x 
is the matrix product of A and x. (1-3) and (1-4) are equivalent 
under the assumption that A . =b, for i =0, 1, ..., k. 

An (n-1)-flat is commonly known as a hyperplane. From (1-4) we see 
that a k-flat is an intersection of n-k independent hyperplanes. If k =n, 
the condition A x = b in (1-4) vanishes. Thus the only n-flat is R” 


itself. A O-flat is just a single real point. 


Example 1-4 In R?, 0,1,2,3-flats are respectively: a single real 


point, a line, a plane, and R? itself. 


A O-flat is a bounded set. A k-flat, where 1 < k < n, is an 


unbounded set. 
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Given k + 1 independent real points in R", there is one and only 


one k-flat containing these real points, and none of smaller dimension. 


1.1.4 Half-flats 


A half k-flat ; in R-. where 1 < k < n, is a k-dimensional unbounded 


convex subset of Ree 
r ={xe€ R | Ax=b and ax < B} (1-5) 


where A x = b defines a k-flat, q is a fixed real n-vector (row) 
which is not in the row space of A, and B is a fixed real number. 


The set 
{x € Re | Ax=b and g x = B} 


is a (k-1)-flat known as the boundary (k-1)-flat of the half k-flat r. 


Example 1-5 In R, half 1,2,3-flats are respectively: a half-line, 
a half-plane, and a half-space. The corresponding boundary 0,1,2-flats 


are respectively: a point, a line, and a plane. 


O 
A half n-flat, 


(x eR" |] a x < B) (1-6) 
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where q is a fixed real n-vector (row), q #0, and f is a fixed real 
number, is commonly known as a half-space. Its boundary (n-1)-flat 
is a boundary hyperplane. An intersection of a half-space (1-6) and 


a half-space, 


{xeR’ | a x28} (1-7) 


is the boundary hyperplane common to both (1-6) and (1-7). 
Let X be a subset of R". A half-space in R" which contains X, and 
whose boundary hyperplane contains one or more real points in X, is 


known as a supporting half-space for X. Its boundary hyperplane 


is known as a supporting hyperplane for X. 


1.1.5 Convex Polyhedra 


n 
A convex polyhedron Pin R is a non-empty intersection of a 


finite number of half-spaces in ae 
n 
P={xeER | Ax < b} 


where A is a fixed real m yxy n matrix and b is a fixed real m-vector 
(column). We shall say that P is an (n,k)-polyhedron if the rank 
of Ais k. 

P may be bounded or unbounded, and may have any dimension between 
0 and n inclusive. Flats and half-flats are special cases of convex 


polyhedra. 
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The structure of a convex polyhedron is best understood by 
considering its faces. A face of a convex polyhedron P is a 
non~empty intersection of P with one or more supporting hyperplanes. 

Let a x < b be any one inequality in the system A x < b defining 


P, The half-space 
tr={xeR > lax <b} 


may or may not be a supporting half-space for P. +r is a supporting 


half-space if 
{x€ P] ax =b} (1-8) 


is non-empty. In this case, (1-8) is a face of P. 


Let S be an arbitrary subset of {1, 2, ..., m}. Let Ag denote 


the submatrix of A comprised of those rows of A whose indices appear 
in S, and let by denote the subvector of b comprised of those components 


of b whose indices appear in S. A subset Pg of P, 


Po ={xeéP | Ag x = bg} 


if non-empty, is a face of P. Furthermore, every face of P may be 
obtained in this manner. 


Lf Py is a face of P, and we let S', where S' ¢ S, identity a maximal 


set of linearly independent rows of Ag, then Pot = Poe Thus we may 


generate all faces P, of P by considering only those sets of indices $ 


S 


which identify linearly independent sets of rows of A. We refer to such 
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sets of indices as face subsets. We include the case S =o , where o 

is the empty set, so that P, which is equal to Ps is a face of itself. 

Since there are only a finite number of ways to form S, P has a 

finite number of faces. The faces of P are themselves convex polyhedra. 
If a face of P has dimension k, we refer to that face as a k-face 

of P. O, I-faces of P, if they exist, are commonly known as vertices 

and edges of P, respectively. 


Tf P. and Por are faces of P and P, ¢ P 


S R gt? then Ps is a subface 


of P Clearly, if S' < S$ then Pg is a subface of P.,. 


st’ S 


A minimal face of P is one which has no proper subface. 
Let P be a convex (n,k)-polyhedron, where 0 < k < n. Hoffman 


and Kruskal [15] show that a face P, of P is minimal iff Ag consists 


S 
of k linearly independent rows of A. Furthermore, they show that if 


Po is a minimal face of P, then 


r = nN _ 


That is, all minimal faces of P are (n-k)-flats. 

The set of all faces of a convex (n,k)-polyhedron P is partially 
ordered by «. The greatest member of this partially ordered set is P 
itself. There is no least member, in general. All faces fall into 
ranks, each rank being composed of all faces of a given dimension. 
The minimal faces of P, which are all (n-k)-flats, constitute the 
lowest rank. The next rank consists of all (n~-k+l)-faces of P, and 
so on. The highest rank consists of P alone, whose dimension lies 


between n-k and n inclusive. 
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Example 1-6 Fig. 1-3 shows graphical representations of several 
2 
convex polyhedra in R and partial ordering diagrams displaying 


the face structure of each. 


The boundary of a face of a convex polyhedron P is the union 
of all proper subfaces of that face. A boundary point of a face 
is a real point in the boundary of that face. A minimal face of 
P has no boundary; that is, its boundary is empty. The boundary 


of P itself is the union of all proper faces of P. 


I~ _- 
P 

x 
~~ 2 27 3 
P 

P P 
P P va ~~ vA aX 
P Ie Paz  Pezzy Peay Pars Porgy Peas 
Fay Pez Paes «xX vi \/ 
Perzz Poisz Peasy = Pavza Pee} 

(a) (b) (c) (d) (e) 


Fig. 1-3 
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If P is a bounded convex polyhedron, then every face of P is 
bounded. Since every minimal face of P is flat, and the only 
bounded flat is a O-flat, it follows that every minimal face 
of P is a O-flat (or vertex), and that P is an (n,n)-polyhedron. 

A bounded convex polyhedron P is equal to the convex hull 
of its vertices (see Hadley [10]). Thus the vertices of P 
completely specify P. 

There is a similar result which applies to any convex polyhedron, 
bounded or unbounded. In order to state it, we first generalize 
the notion of a minimal face. 

A face of a convex polyhedron P may or may not be equal to the 
convex hull of its boundary. A face of P is said to be reducible 
if it is equal to the convex hull of its boundary, and irreducible 
if it is not. Flats and half-flats are irreducible faces. 

Klee [19] shows that a face of a convex polyhedron P is 
irreducible iff that face is a flat or a half-flat. He shows 
further that P is equal to the convex hull of its irreducible faces. 

A minimal face of P is irreducible, but an irreducible face 
of P is not necessarily minimal. The converse does not hold 
because a half-flat is irreducible but not minimal. The notions 
of an irreducible face and a minimal face coincide in the case of 
bounded convex polyhedra since no face of a bounded convex polyhedron 
can be an (unbounded) half-flat. 

If an irreducible face of P is a half-flat, then its boundary 


is a flat which is also an irreducible face of P. A maximal irreducible face 
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of P is one which is not contained in another irreducible face of P. 
Klee [19] proves the following sharper result: 
A convex polyhedron P is equal to the convex hull of its maximal 


irreducible faces. 


Example 1-7 The minimal, irreducible, and maximal irreducible faces 


of the convex polyhedra represented in Fig. 1-3 are tabulated below: 


Maximal 

Fig. | Bounded/ Minimal Irreducible 
Irreducible 
1-3 Unbounded Faces Faces 
Faces 
(a) Bounded P 
(b) Unbounded P 
P P P Pp P Pp 

a a fy? *¢23 (13? Peay cay? Peay 


(d) Bounded Pras2}? Praia}? Prai2}? Peai3}? Fra,2}? *(h3)? 


P23} eee Fo 3} 


(e) Unbound ed Praia}? Pro 3} Pray Pray}? Pray? Pe 3} 


Pea ay? *¢2,3) 
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A class of convex polyhedra of special interest are those 
arising in connection with linear programming problems. A linear 


programming problem has the form: 


Maximize c x 
Subject to Ax < b (1-9) 
x 2 0 (1-10) 


where c is a fixed real n-vector (row), A is a fixed real m y n matrix, 
and b is a fixed real m-vector. The linear function cx is known as 

an objective function. A linear programming problem may instead 
involve minimizing an objective function. By simply negating c we 

may change such a minimization problem into a maximization problem. 


The set P of real points satisfying the inequalities (1-9) and (1-10), 
P=(xeR |Ax<b and x 20} (1-11) 


is a convex polyhedron. We shall refer to P as a linear programming 
polyhedron. The condition x = 0 gives P some special characteristics 
which we now examine. 

Suppose the inequalities (1-9) and (1-10) are consolidated and 


(1-11) is rewritten as 


P={xe R7 | Alx< b'} 


23 


where A' contains as a submatrix the n y n identity matrix (negated). 
The rank of A' is n, and it follows that P is an (n,n)-polyhedron, 
all of whose minimal faces are O-flats. 
What about the irreducible faces of P? It can be shown that 
{x € a | x > 0}, known as the nonsnegative orthant of Re contains 
no l-flat. It follows that the non-negative orthant of R" contains 
no k-flat, for 1<k <n, and no half k-flat, for 2<k <n. Since 
P is wholly contained in the non-negative orthant of ae P contains 
no k-flat, for 1 <« k <n, and no half k-flat for 2 <¢k <n. Thus 
the irreducible faces of P, necessarily being flats or half-flats, 
must be O-flats or half 1-flats, 
P is equal to the convex hull of a finite number of O-flats and/or 


half l-flats. 


L236 Integral Convex Sets 


Let J be the ring of integers. Let J” be the n-dimensional 
module over J. We refer to x ¢ a6 as an integer point or an 
integer n-vector. 

For any subset Xc R we define I(X) to be the set of integer 


points contained in X, 


L(X) = {xe J” | x © X}. 
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A convex subset X ¢ R is said to be integral (or to have the 
integral property) iff X is the convex hull of the set of integer 
points contained in X, that is, iff X = H(I(X)). Since I(X) < X 
it follows that H(I(X)) ¢ H(X) = X. Thus X is integral iff 
X ¢c HCI(X)), that is, iff every x € X is a convex combination of 
integer points in X. 

We are primarily interested in convex polyhedra having the 
integral property or, more simply, integral convex polyhedra. 

In Chapter 2 we develop necessary and sufficient conditions for 
a convex polyhedron to have the integral property. 
As the following example shows, not every integral convex 


set is an integral convex polyhedron. 


0 
Example 1-8 Let X ={x , Pa ae --- } be the infinite subset 


of i defined inductively as follows: 


(0, 0) 


” 
Mt 


Cae x + (n+1, 1). 


Fig. 1-4 shows the first few members of X and a portion of the boundary 

of H(X). H(X) is an integral convex set because every x ¢€ H(X) is a convex 
combination of integer points in H(X). Clearly H(X) is equal to the 
intersection of an infinite number of half-planes and therefore is not 


a convex polyhedron. 
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Fig. 1-4 


If X is a convex subset of R" then there exists a unique convex 
subset X' of R” defined by X' = H(I(X)). In words, X' is the convex 
hull of the set of integer points contained in X. We point out two 
extreme cases. First, if X is an integral convex set then 


x! 


H(I(X)) =X. Second, if X contains no integer points then 


x! 


H(~) =o, where ~ is the empty set. 


The following lemma is a direct consequence of our definitions. 


Lemma 1-1 


Let X and X' be convex subsets of R”. Then xX! = H(I(X)) iff 


(i) X' is integral 
and 


(ii) IK") = I(x). 


27 


Example 1-9 Consider the following system of linear inequalities: 


(a) 4x < 15 
(b) 2x < 9 (1-12) 


(c) -3x < cl 


The shaded interval of Fig. 1-5(i) represents the set of real 


points x which satisfy (1-12). 


(i) (ii) 


Fig. 1-5 


Some of these real points are integer points and these are indicated 
with dots. In Fig. 1-5(ii) the shaded interval is the convex hull 
of these integer points. It represents the set of real points 


which satisfy the following system of linear inequalities: 


28 
(a) x < Le | =3 


(b) -x < |= 


where | J denotes the greatest integer less than or equal to 


-1 


the real number A. The shaded interval of Fig. 1-5(i) is an 

example of a convex polyhedron. The shaded interval of Fig. 1-5(ii) 

is an example of an integral convex polyhedron. The process of 
transforming the former into the latter is an example of integralization. 


O 


Example 1-10 Consider the following system of linear inequalities: 


(a) 2x, + 2x5 < 15 


1 
(1-13) 


(b) “4x, - 4x5 < 718 


The shaded region of Fig. 1-6(i) represents the set of real points 

x = (x, > x5) which satisfy (1-13). It contains a countably infinite 
mumber of integer points, some of which are indicated with dots. In 

Fig. 1-6(ii) the shaded region is the convex hull of these integer points. 
It represents the set of real points which satisfy the following system 


of linear inequalities: 
15) _ 
(a) gr ie el Ger ia 


@) xo xy < [AB] = 


In Fig. 1-6(i) we have an example of an unbounded convex polyhedron, 
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and in Fig. 1-6(ii) an example of an unbounded integral convex polyhedron. 


Again, the transformation of the former into the latter is an example 


of integralization. 


+ 


i ; 
See ene eee eee 


(ii) 


(i) 


ig. 1-6 


F 


Consider the system 


Example 1-11 


+ Ty < 29 


(a) x) 


(1-14) 


(b) AX, = 5X5 < 12 


“7%, - 4X, < 727 


(c) 


30 


Fig. 1-7(i) illustrates the convex polyhedron defined by (1-14), 
while Fig. 1-17(ii) illustrates the integral convex polyhedron 


which is the convex hull of integer points contained in the former. 


(i) (it) 


Fig. 1-7 


The integral convex polyhedron of Fig. i1-7(ii) is defined by the 


system of inequalities, 


WA 
ies) 


(a) Xy 
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If we were to imagine the x plane as a peg~-board, with pegs 


3 
placed in all integer point positions, then a rubber band in tension, 
made to conform with the boundary of the shaded region of Fig. 1-7(i), 
would, when released, assume the form of the boundary of the shaded 
region of Fig. 1-7(ii), assuming that the rubber band remains in 
tension at the completion of its contraction. In its final state, 
the rubber band is supported entirely by pegs. Every vertex of 
the region it encloses is an integer point. This is a physical 
method for integralizing the bounded convex polyhedron of this 
example. 

In three dimensional space, a rubber sphere made to conform 
with the boundary of a bounded convex polyhedron and enclosing 
small marbles somehow fixed in integer point positions would do 
the job. In higher dimensional space it is difficult to imagine 
any physical process for accomplishing integralization. What is 
nedded, of course, is a mathematical procedure for performing 


the integralization transformation. oO 


In the chapters to come we shall be primarily concerned with 
two questions: 

(i) Given a convex polyhedron, how does one tell if it is an 
integral polyhedron? 

(ii) Given a non-integral convex polyhedron, how does one 


integralize it? 
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We answer the first question in Chapter 2 and we answer the 


second question for a restricted class of polyhedra in Chapters 


4, 5 and 6. 
1.3 Integer Linear Programming and Related Research 


In this section we seek to accomplish two goals. Our first 
goal is to establish the relevance of integral convex polyhedra 
and the process of integralization to integer linear programming. 
Our second goal is to describe some of the research done by others 


concerning integral convex polyhedra and integralization. 


1.3.1 Integer Linear Programming 


As we stated in the preceding section, a linear programming 


problem has the following form: 


Maximize c x 


Subject to Ax <b (1-15) 


x 2 O 
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where c is a fixed real n-vector, A is a fixed real m y n matrix, 
and b is a fixed real m-vector. The inequalities in (1-15) define 
a convex polyhedron P. 

A commonly used algorithm for solving (1-15) is the simplex 
method developed by Dantzig [5]. The fundamental idea of the 
method is: if an optimal solution to (1-15) exists then some 
vertex of P is an optimal solution. Beginning with any vertex 
of P, the simplex method generates a sequence of vertices of P 
for which the objective function cx takes on non-decreasing values. 
The algorithm terminates when an optimum vertex is found, or 
gives an indication that no optimal solution exists. 

The simplex method has proven to be an efficient algorithm 
and has found many applications in operations research, management 
science and economics. 


An integer linear programming problem has the form: 


Maximize cx 


Subject to Ax <s b 
(1-16) 
x 2 O 


x an integer vector 


where c is a fixed integer n-vector, A is a fixed m yy n matrix, 
and b is a fixed integer m-vector. There are variations to (1-16), 


notably the more general formulation in which all constants are 
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allowed to be real. In another variation, only a fixed subset of 
the components of x are required to be integers; the remaining 
components are allowed to be real. We shall take (1-16) to be 
the standard integer linear programming problem. 

Again, the inequality constraints of (1-16) define a convex 
polyhedron P. The simplex method, when applied to (1-16), does 
not work in general simply because it is no longer true that if 
an optimal solution to (1-16) exists then some vertex of P is an 
optimal solution. Some or all of the vertices of P may not be 
integer points. Said in another way, optimal solutions to (1-16), 
if they exist, may be embedded in the interior of P and not 
accessible to any vertex generating algorithm. 

If P happens to be an integral convex polyhedron then, as we 
shall see in Chapter 2, every vertex of Pis an integer point and 
the simplex method will produce an optimal solution to (1-16) 
for any objective function cx. 

Suppose P is not an integral convex polyhedron, as is usually 
the case. If we integralize P to obtain the integral convex 
polyhedron P', where P' = H(I(P)), then the simplex method applied to 
P' will produce an optimum solution to (1-16) for any objective 
function cx. We acknowledge the fact that there is a certain amount 
of naivety in this approach to integer linear programming. This 
stems from the fact that P' ts potentially a very complex object — 


complex in terms of the number of inequalities needed to define Pi — 
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and it certainly seems inefficient to construct P' in order to 
discover just one of its vertices. 
However, we offer the following observation, made by Gomory, 
in support of the idea of integralization as an important subject of 
research. We quote from Gomory [11]: 
"Since any algorithm for the integer programming problem, whether 
related to linear programming, branch and bound, exhaustive search, 
or whatever, must end up finding a vertex of P', information on 
P' seems relevant to any approach to the integer programming problem. 
Yet information about P' is very difficult to obtain." 
Although the last sentence of this quotation does not contribute 
to our argument, we include it to avoid mis-representation by omission. 
It has been our goal to gain some insight into the structural 
relationship between P' and P by focusing on the integral property 
and the process of integralization for convex polyhedra in general, 
not just linear programming polyhedra. We believe that such insight 
can contribute to the improvement of existing integer linear 


programming algorithms and the development of new ones. 


1.3.2 Related Research 


In this subsection we describe some of the known results related 


to integral convex polyhedra and integralization. 
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The major results concerning integral convex polyhedra known 
to us are contained in the paper by Hoffman and Kruskal [15]. 
They consider convex polyhedra defined by system of linear inequalities 
Ax < b in which all constants are integers, and define such a convex 
polyhedron P to be integral iff every minimal face of P contains 
an integer point. If the minimal faces of P are vertices then 
this definition says that P is integral iff every vertex of P 
is an integer point. As we shall see in Chapter 2, their definition 
is equivalent to our own. 

Hoffman and Kruskal obtain the following result concerning 
the integral property for convex polyhedra. Let A be an integer 
my nm matrix of rank k. Let S denote a subset of row indices 
of A and let Ag denote the submatrix comprised of rows of A identified 
by S. Let ged [Ag] denote the greatest common divisor of all 


|S|-order minors of A Then Ax < b defines an integral convex 


Ss? 
polyhedron P(b) for every integer m-vector b iff 


ged [A,] =1 (1-17) 


for all k-subsets S of linearly independent rows of A. We observe 
that for a particular integer m-vector b, condition (1-17) is 
sufficient for P(b) to be integral but not necessary. 

Hoffman and Kruskal also obtain a result concerning the integral 


property for linear programming polyhedra. An integer m y n matrix 
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is said to have the unimodular property iff every minor (of every 
order) of A has value 1, 0, or ~1. Entries in a unimodular matrix 
are therefore restricted to be 1, 0, or -1. They show that the 
inequalities Ax < b and x 2 O define an integral convex polyhedron 
Q(b) for every integer m-vector b iff A is a unimodular matrix. 
We observe that for a particular integer m-vector b, the condition 
that A be unimodular is sufficient for Q(b) to be integral but 
not necessary. Unimodular matrices occur in the integer linear 
programming formulations of several combinatorial optimization 
problems, including the network flow problem and the maximal 
matching problem for bipartite graphs. 

Concerning integralization, Gomory in [10] are [11] gives 
a method for integralizing a corner polyhedron P — that is, 
a polyhedron defined by Ax < b where A is an integer n xy n matrix 
of rank n and b is an integer n-vector. Basically the method 
involves generating a large but finite number k Of auxiliary 


inequalities, 


cx 2 1, i#1, 2, ...,k 


. 


i : ; 
where c , for i=l, 2, ..., k, is obtained as an integer solution 
to an equation over a finite Abelian group derived from the matrix A. 


Gomory shows that these inequalities define a polyhedron whose vertices 
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correspond one-to-one to inequalities defining P', where P' = H(I(P)), 
and indicates how to obtain the latter from the former. Thus the 
problem of integralizing a corner polyhedron becomes one of 
finding all vertices of an auxiliary polyhedron, which can be 
accomplished by linear programming methods. 

We mentioned previously that the end product P' of integralizing 
a convex polyhedron P is potentially a very complex object. Rubin [22] 
obtains an interesting quantitative result concerning this complexity. 
Let Ays Ags Ags sees aps vee be the infinite sequence of integers 


defined by the recurrence, 


i 


ay ae + a2 for k 23 


=l, ay = ll. 


"4 
This is the well known Fibonacci sequence. Rubin defines an infinite 
sequence Pi» Po» P3> eee, Poe «e- of bounded (2,2)-polyhedra based 


on the Fibonacci sequence as follows. Phe for k 21, is defined by 


the three inequalities, 
+ x 
Fa *1 Fatt *2 SF 2k+1 
x, 2 0, Xo 2 0. 


Pie for k 2 1, has 3 vertices and 3 edges. He shows that Pu where 


Pe = H(I(P, )), for k 2 1, has exactly k+3 vertices and k+3 edges. Thus 


39 


for a 2-dimensional bounded polyhedron P, it may require 
arbitrarily many inequalities to define P' = H(I(P)), even 


though P is detined by just three inequalities. 
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CHAPTER 2 


THE INTEGRAL PROPERTY IN CONVEX POLYHEDRA 


In this chapter we develop several results concerning convex 
polyhedra having the integral property. Using our definition of 
an integral convex set X as one which is the convex hull of the 
set I(X) of integer points in X, we derive necessary and sufficient 
conditions for convex polyhedra to have the integral property. 

The chapter is divided into three sections. [In the first 
we deal with integral flats, in the second with integral half- 


flats, and in the third with integral convex polyhedra in general. 


2.1 Integral Flats 


A k-flat 7 in R" is generated by taking all affine combinations 
of an independent set of k+l real points. Lemma 2-1 tells us when 


7 is integral in terms of its generating set. 


Lemma 2-1 
A k-flat , is integral iff mn is generated by an independent set 


0 
{x , oe Setar nS of integer points. 
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Proof 
We first prove necessity of the condition. Assume 7 is an integral 
k-flat. Since 7 is a k-flat, 7 is generated by an independent set Y of 
k+l real points. Since 7 is integral, each real point in Y is a 
convex combination of integer points in I(r). The dimension of I(r) 
is at most k. If I(t) has dimension less than k, then Y has dimension 
less than k, which is a contradiction. Thus I(t) has dimension k and 
m contains an independent set of k+l integer points which generates r. 
In order to prove sufficiency of the condition, assume 7 is 


a sieves <) of integer points. 


; 0 
generated by an independent set {x , x 
We shall show that every x ¢ mn is a convex combination of integer 
points in y. It will then follow that nm is integral. 


Let x be an arbitrary real point in mn, 


0 1 k 
x = Xo x + Ny xX + eee + oF x 


where No» Ni» Sereks Me are real numbers and 


Ag tA, tote tA ed. 


Equivalently we may write 
0 2 
x =x + h, (xix?) + No (x ey totee ft hy (Sox®). (2-1) 


Let functions ho: R»J and h R » J be defined as follows. 


1° 


Let ho @,) be the greatest integer less than or equal to AS and 
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let h, @,) be the smallest integer greater than hie (Thus 

ak 
h, @,) - hoQ,) = 1, for all real Ae) let {0,1}” be the set 
of all binary k-tuples, a representative member of which is 
denoted by b = (b> by» wens bs We claim that x is a 


. : k, : 
convex combination of the 2° integer points, 


k 
Ze = Ze + y. h, O,;) cha) all be (0,1}* 
i 


i=1 


in ~. (Geometrically, these ok integer points are vertices of a 
k-dimensional cube in 7 and we claim that x lies in this cube.) 

In order to exhibit this convex combination we define functions 
Bq? R»R and 84: R + R as follows: 


8, A,) =A; ~ho@,)- 


We then define function &,° a 4 R as_ follows: 
k 


BO) = al &, O,)- 


i=l 


We claim that 


(2-2) 
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where 
k 
O<g,Q) <1 for all b € (0,1) 
and 
y gQ) =1. 
befor)" 


In (2-2), % = O,> hoe ees ri is the real k-tuple of coefficients 
appearing in (2-1). 

In order to prove this claim we first observe that (2-3) holds 
since 0< By A,) <1 and O< 8, A,) <1 for all real numbers hy 


Next we prove (2-4) by induction on k. For k = 1 we have 
Bo) + 8,0.) = h,A,)-hy@,) = 1 


Assume (2-4) is true for some k, 1 <k <n. Then 
k+l 


: [ ll ee ] 


pre (0,4) 


k 
Exuere - 2p | "pee La [ il m0] 


bo Ore) + 8] Oy? 


BOD “ho Aap = 1 


(2-3) 


(2-4) 
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Finally we prove that (2-2) is valid. We have 


ie 

=! y oR 0 v o _ iL 0 

-/ os g,) | ge ‘ 850 i dh A) Grex )| 
b ¢ {0,1} be 09 a 


yy 


k 
= 30 my y r x BO) hy | Cte 


The expression within brackets may be expanded as 


} Bp ~) A : A, ) 


be fo.1y" 
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k 
=[ a ,) boA,) + 8,0) 1, O,)]- FD [ TL Bp, a) | 
| b € (0,1)* 

j#i 


By, Ho®,) + 8, 0,) hLO,) 


(hy Aj I-Az) AoA.) + A,-ho®,)) hLA,) 


r, (hy Q,) ~ ho®,)) 


=-h, 
L 


Thus i: 


ee i r, Gtx?) 
isl 


y 8, O) x? 
e (0,1}* 


o 


and we have expressed an arbitrary x € y as a convex combination of 

integer points in tr. It follows that 7 is integral. oO 
Let A be an integer (n=-k) y m matrix of rank n-k and let b 

be an integer n-vector. Our next objective is to determine what 

conditions must be imposed on A and b in order for Ax = b to 

define an integral k-flat. The preceding lemma allows us to 

rephrase the question as: What conditions must be imposed 

on A and b in order for the solution space of Ax = b to contain 


an independent set of k+l integer points. 
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We first introduce the following notation. If B is an integer 
ry s matrix of rank r, where 1 < r < s, then we define gcd [B] 
to be the greatest common divisor of all ) ro} _oyder minors 
of B, where @) is a binomial coefficient. Notice that if r =s 
then B is a square matrix and gcd [B] = |det B], while if r =1 
then B is a row vector and gcd [B] is the greatest common divisor 


of all components of B. 


Theorem 2-1 
Let A be an integer (n-k) y n matrix of rank n-k and let b 


be an integer (n-k)-vector. Then 


> 

bad 
It 
o 


defines an integral k-flat ry iff 


gcd [A] = gcd [A ib]. 


Proof 

Assume that A x= b defines an integral k-flat 7. By Lemma 2-1, 
7 is generated by an independent set of k+l integer points. Let 
a be an integer point in 7. Then Ax? =b. Let S denote an 
arbitrary (n-k-1)-subset of column indices of A and let ae be the 
(n-k) yx (n-k-1) submatrix of A comprised of those columns of A 


identified by S. Since 


47 


det fan b] = det ta® oar) 


S. t 
it follows that det [A !b] is an integer combination of (n-k) Be stdar 


minors of A. Thus 


St Bl 


ged [A] | det [A 
Since this holds for any S, we conclude that 


gcd [A] = ged [Ai b]. 


Now we assume that 


gcd [A] = gcd [A jb] 


and prove that Ax =b defines an integral k-flat rn. 
According to Jacobson [18], there exists an integer (n-k) x (n-k) 
matrix R with det R = +1 and an integer n xy n matrix C with det C= + 1 


such that 


A=RAC 


where A is an integer (n-k) yx n diagonal matrix, 


° 
eee ee 
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Furthermore, Jacobson [18] proves that 


ged [A] = ged [A]. 
Tt is easily seen that 
ged [A] = ay ay ttt ale 
We may write A x =) as 
RACx=b 
or as 
A y=b 
where 
y=cCcx 
ee ae 


Tf we let C' be the following integer (n+1) y (n+l) matrix 


\ 

lo 
. ' 
Cc’ = I 

—— ee ee tO 

0 il 

i 
we see that det C' = det C = #1 and 


[Aib] =R [A‘tb] cl. 
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Again, it follows from Jacobson [18] that 


ged 


We now have 


fA: 


vee 0 
a 
eed ‘ : 
6 2 
Foek 
Thus 
aA 
ay a, 
a 38, 
and 


| 0 = ced 


<; ate | ay BD ath 
a" ak | a ay ttt b 
a, | >, 
a, | b, 
: : 
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The k-flat defined by A y= b contains the following independent 


set of k+l integer points, 


b,/a, by /a, bi fay b,/a, 
by/ay by/ay bo/a, by/ay 
bk! nHk ake ak bk Aa-k oes bak Fank 
0 1 0 0 
0 0 1 0 
0 0 0 1 


Therefore a is generated by these k+l integer points and, by 
Lemma 2-1, 7 is integral. Since the k-flat 1 defined by A x =b 
is the image of - under the integer transformation gt. mr is 
also an integral k-flat, which is what we set out to prove. 0 
Corollary 2-1 

A x =b defines an integral k-flat iff A x = b has an 


integer solution. 
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Proof 
The necessity of the condition follows from Lemma 2-1. 
To prove sufficiency, assume A x = b has an integer solution 


0 
x. Then, as in the necessity part of the proof of Theorem 2-1, 


gcd [A] = gcd [A ib]. 


Thus, by Theorem 2-1, Ax 


b defines an integral k~flat. 
O 

Suppose that we have a system A x = b of n-k linear equations 
which defines an integral k-flat. The following corollary shows 
that any (n-k')-subset, k < k' < n, of these linear equations 
defines an integral k'-flat, Let S be an arbitrary (n-k')~subset 
of {1,2,..., n-k}, where ke k' <n. Let A, denote the submatrix 
of A comprised of those rows of A whose indices appear in S and 
let be denote the subvector of b comprised of those components 


of b whose indices appear in S. 


Corollary 2-2 


If A x = b defines an integral k-flat then Ag x = be defines 


an integral k'-flat. 


Proof 


Suppose A x = b defines an integral k-flat. By Corollary 2-1, 


i . 0 0. ; 
A x = b has an integer solution x . But x is also an integer 
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solution to Ag xX =D Thus, again by Corollary 2-1, A, x =b 


s* s S 


defines an integral k'-flat. 


2.2 Integral Half-flats 


We now consider the conditions under which a half-flat is 


integral. Let 7 be a half k-flat in R”, 
yr = {xe R® | Ax =b and ax < B} 


and let m be its boundary (k-1)-flat. Again, we consider all 


constants to be integers. 


Theorem 2-2 


A half k-flat - is integral iff its boundary (k-1)-flat 7 is 


integral. 


Proof 


Assume 7 is integral. Since -t is not empty, I(r) is not empty. 
If m7 contains an integer point then, by Corollary 2-1, , is integral. 


Otherwise 


a x< B for all x € I(r) 
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or 


a x < Berl for all x € I(r). 


Thus A x = b and g x < Brl define a new half k«flat +' such that 
r' cere and I(r¢') = I(,). But this contradicts our assumption 
that - is integral. Thus 7, contains an integer point and 7 is 
integral. 

Now assume that 7 is integral. We shall argue that 7 is 
the convex hull of a countably infinite number of integral 
(k-1)-flats, and thus 7 is itself integral. 


Let 


It can be shown that 7 is the convex hull of the set of (k-1)-flats 
{n> 2 TTo2 eee}. We shall show that each (k-1)-flat 7; in the 


set is integral. 
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First of all, 709 is integral by assumption. We then have from 


Theorem 2-1 that 


ged = gcd 


| 

A A |b 
| 
| B 


th 
We next observe that each (n-k+l)” <order minor of 


which includes the right hand column differs from the corresponding 


minor of 


by a multiple of A. It follows that 
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and by Theorem 2-1, 7; is integral. Then, since each x € y is a 
convex combination of real points in To» Ty Ts eee 5 @ach of 
which is a convex combination of integer points in 7, it follows 
that each x € 7 is a convex combination of integer points in +. 


Thus , is integral. 


2.3 Integral Convex Polyhedra 


We are now prepared to derive some conditions under which a 
convex polyhedra has the integral property. Let P be a convex 


polyhedron in Rey 
n 
P={x¢eR | Ax <b}. 


As before, we assume all elements of A and b to be integers. 


Theorem 2=3 


P is integral iff every face of P is integral. 


Proof 
Suppose every face of P is integral. Since P is the convex hull 


of its irreducible faces all of which are integral, P is itself integral. 
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Conversely, assume P is integral. If P is a flat then its 
only face is P itself which is integral by assumption. Otherwise 


P has one or more proper faces. Let Py be any proper face of P, 


Pg ={xcP | Ag x = bg} 


where S specifies a subset of the inequalities defining P. 


Let x € Poe Since x ¢ P, x is a convex combination of integer 


2 
points x! KX 5 cee, x" e P, 


We claim that x, ne eee, x" € Pos that is 


A, x =b for i =1,2,..., m. 
Otherwise 


i ‘ 
Ag x < by for some i 


and 
m m 
= i_ 1 
Ag % = Ag ye = Vy, cy Ag ® < bg 
i=l 


which contradicts the fact that x ¢€ Poe Thus x is a convex combination 


of integer points in Po? and P, is integral. 


° O 
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The following corollary gives a necessary and sufficient 
condition on the system of inequalities A x < b defining P for 


P to be integral. 


Corollary 2-3 


P is integral iff 


ged [A] = ged [Ag } bg] 


s 


for all face subsets S. 


Proof 


Assume P is integral. By Theorem 2-3 every face Py of Pis 
integral. Let Tg be the flat defined by Ag x = bes Since Po Sg» 


Tg must contain an integer point, and by Corollary 2-1 and Theorem 2-1, 


gcd [Ag] ged [A, ‘bol. 


S 


Now assume that 


gcd [Ag] = gcd [A, i bo] 


for all face subsets S of P. Consider the irreducible faces of P. 


An irreducible face P, of Pis either a flat or a half-flat. If Po 


is a flat then P. is integral by Theorem 2-1. If Po is a half-flat 


then its boundary flat is an irreducible face of P which is integral 
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by Theorem 2-1. Then, by Theorem 2-2, P, is integral. Since P 


S 
is the convex hull of its irreducible faces, all of which are 


integral, P is integral. oO 


We may further restrict the condition of the preceding corollary. 


This is done in the following. 


Corollary 2-4 


P is integral iff 


ged [A,] = ged [A ‘bol 


S S 


for all minimal-face subsets S. 


Proof 


Necessity follows directly from Corollary 2-3. In order to prove 


that the condition is sufficient, assume that 
ged [A,] = ged [A, ib] 


for all minimal-face subsets S. Consider the irreducible faces of P. 


An irreducible face Py of P is either: (i) a flat which is a minimal 


face of P, or (ii) a half-flat whose boundary flat is a minimal face 


of P, In the first case, P, is integral by Theorem 2-1. In the second 


) 
case, the boundary flat of Po is integral by Theorem 2-1, and thus Po 
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is integral by Theorem 2-2, Since P is the convex hull of its 


irreducible faces, all of which are integral, P is integral. 


Corollary 2-4 tells us that P is integral if all of its 
minimal faces are integral. If P is an (n,k)-polyhedron then 
the rank of Ain Ax <b is k and the minimal faces of P are all 
(n-k)-flats. If k =n then the minimal faces of P are all 
vertices. For example, this is the case for all linear programming 
polyhedra since the non-negative condition x = 0 insures a rank n 
coefficient matrix. In such cases, Corollary 2-4 asserts that P 


is integral iff every vertex of P is an integer point. 
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CHAPTER 3 


THE TOTALLY INTEGRAL PROPERTY FOR CONVEX POLYHEDRA 


In this chapter we define a new property which we call the totally 
integral property. We show that the class of totally integral polyhedra 
is a proper subclass of the class of integral polyhedra. We then 
consider the class of linear programming polyhedra having the totally 
integral property and show that this class can be characterized 
algebraically using some of the group theoretic ideas developed by 


Gomory [9]. 


3.1 Totally Integral Convex Polyhedra 


We begin by defining the totally integral property. Let P be 
a convex polyhedron defined by the system of inequalities A x < b, 
where, as usual, all elements of A and b are integers. Let S be 
a subset of row indices identifying a linearly independent set 


of rows of A. Then the system of equations A, x = b, defines a flat 


S 
which we denote by TIge We shall say that Tg is a flat of P and that 


S is a flat subset Notice that a face subset is always a flat subset, 
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but that a flat subset need not be a face subset. It Tg is a flat 


of P then it is possible that Pos where 
Po = Pang 


is empty and thus not a face of P. 
A convex polyhedron P is said to be totally integral iff 


every flat Tc of Pis integral. 


Example 3-1 Fig. 3-1(a) illustrates a totally integral convex 
polyhedron. Fig. 3-1(b) illustrates an integral convex polyhedron 


which is not totally integral. oO 


(a) (b) 


Fig. 3-1 
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The following lemma gives us a necessary and sufficient 
condition for P to be totally integral, expressed in terms of 


A and b. 


Lemma 3-1 


P is totally integral iff 


ged [Ag] = gcd [A, jb,] 


for all flat subsets S. 


Proof 


The lemma follows directly from the definition of the totally 


integral property and Theorem 2-1. 


O 


Suppose matrix A has rank n-k, where 0 < k en. We shall say that 
a flat Tg of P is minimal iff S identifies n-k linearly independent 
rows of A. Thus a minimal flat of P is a k-flat. We may restrict 


the condition of the preceding lemma as follows. 
Lemma_ 3-2 
P is totally integral iff 
ged [Ag] = ged [Ag ib] 


for all minimal-flat subsets S. 
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Proof 


Necessity follows directly from Lemma 3-1. To prove sufficiency, 


assume 


gcd [A,] = ged [A, ib] 


for all minimal-flat subsets S. Let R be any flat subset. We may 
augment R with additional rows of A, if necessary, to obtain a 


minimai-flat subset S such that Ro S. By assumption, 


ged [Ag] = gcd [Ag }bg] 


and by Corollary 2-2 and Theorem 2-1, 
= : 7 
gcd [A] gcd [A, 1b, ] 


It follows from Lemma 3-1 that P is totally integral. oO 
We see immediately that a totally integral convex polyhedron P 
is always integral. Every minimal face of P is a minimal flat of 
P. If P is totally integral then every minimal face of P is integral 
and, by Corollary 2-4, P is integral. 
The converse is not true ™ an integral convex polyhedron P 
need not be totally integral. It is possible that every minimal 
face of P is integral, but that one or more minimal flats of P 


which are not faces do not have the integral property. 
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Suppose we have an integer m y n matrix A and we let T(A) be 


the set of integer m-vectors b such that A x < b defines a totally 


integral convex polyhderon. We show that T(A) forms a group under 


vector addition. 


We first observe that T(A) always contains the all zero m-vector. 


Applying Lemma 3-1 with b = 0, since 


gcd [Ag] = ged [A, +0] 


for all flat subsets S, it follows that 0 ¢ T(A). Furthermore, T(A) 


rae 1 2 : 
is closed under vector addition. For if b and b are two integer 


m-vectors in T(A) then 


and 


a the 
gcd [Ag] = gcd [Ag ibe] 
gcd [A,] = ged [A He] 

S s''s 


for all flat subsets S, and it follows that 


_ oe ee 
gcd [Ag] = gcd [A, ibe + bg] 


2 , 
for all flat subsets S. Thus be +b ¢€ T(A). Moreover, T(A) contains 


inverses for each of its members. If b ¢€ T(A) then 


and 


gcd [Ag] = 


I 
gQ 
Q 
[ary 
= 
> 
n 
o 
he 


gcd [Ag] = 


I 
ag 
Q 
jon 
— 
> 


for all flat subsets S. Thus <b ¢€ T(A). 
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In the following section we completely characterize the group 


T(A) for linear programming polyhedra. 


3.2 Totally Integral Linear Programming Polyhedra 


We now focus our attention on those convex polyhedra contained in 
the non-negative orthant of RY — that is, convex polyhedra P defined 
by 

Ax<b 

(352) 
x >» 0 

where again A is an integer m y n matrix and b is an integer m-vector. 
Given matrix A, let T(A) denote the set of integer m=vectors b for 
which (3-1) defines a totally integral convex polyhedron. In the 
following, we obtain a complete characterization of T(A), using 
group theoretic methods originally developed by Gomory [9]. 
Specifically, we show that T(A) may be characterized implicitly as 
the set of integer solutions to a single equation expressed over 


a finite Abelian group. 
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Writing the constraints (3-1) defining P as 
A b 
jn-=—-| X SS ae 


where I is ann y n identity matrix, we see that the coefficient matrix 


B has rank n. It follows that all minimal flats of P are O-flats 


-I 


(or points), as are all minimal faces of P. From Lemma 3-2, P is 
totally integral iff all minimal flats of P are integer points. 
Also, from Corollary 2-4, P is integral iff all minimal faces of P 
are integer points. 


We reformulate the constraints (3-1) defining P in the following 


[Ai T] B =b (3-2) 
Zz 


x20, z20 (3-3) 


equivalent form, 


where I is an my m identity matrix and z is an m-vector of slack 
variables. A basic solution to (3-2) is one which is described 
as follows. Let B be any basis of [A! I], i.e. a non-singular 


m y m submatrix of [A: I], and let Xp be the corresponding 
m-subvector of J . Let N be the remaining m y n submatrix 


of [AiI] and let xy be the corresponding n-subvector of Bl . 
Zz 


We may write (3-2) as, 
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The unique solution to (3-4) given by 


Xy = 0 


is a basic solution to (3-2). If x, 20 then the basic solution 


= B 


aps is a basic feasible solution to (3-2), and B is a feasible 


basis for (3-2). Notice that whether or not a basis B of [A:I] is 


feasible depends upon the particular right hand side vector b in 


(3-2). 


It is easily seen that O-flats of P correspond to basic solutions 
of (3-2) andthat O-faces (or vertices) of P correspond to basic 
feasible solutions to (3-2). Thus, P is totally integral iff every 
basic solution to (3-2) is an integer point, and P is integral iff 
every basic feasible solution to (3-2) is an integer point. 

We now take a slightly different point of view in expressing 
the condition for P to be totally integral. Let B be a basis 


and let M(B) denote the set of all integer combinations of columns 


x 


of B. Then the basic solution oe is an integer point iff the 


0 


(3-4) 
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solution * to Bx, = b is an integer point or, equivalently, iff 

ntm ' : 
b € M(B). There are at most ( as ) bases of [A!I], which 
we denote by Bi Bos se Bye It follows that P is totally 


integral iff 
be M(B, ) for te 12 y seas 
or, equivalently, iff 


k 
a a M(B, )- 


i=l 


Recalling our defining of T(A) as the set of integer m-vectors 
b for which P — defined by (3-2) and (3-3) with A fixed — 


is totally integral, we see that 


T(A) = a M(B, ). 
i=l 
We note that M(B,), for i = 1,2, ..., k, forms an Abelian group 
under vector addition. Moreover, this group is a subgroup of 
the Abelian group M(I) under vector addition, where I is the 
m x m identity matrix. It follows that T(A), as an intersection 
of subgroups of M(I), is also a subgroup of M(I). We summarize 


these results in the following. 
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Lemma 3-3 


Given a fixed integer m y n matrix A, the set T(A) of all 


integer m-vectors b for which the constraints 


[Ait] |=] =b 
VA 


x 2 0, z>2O 


define a totally integral convex polyhedron is a subgroup of M(I). 


Moreover 


T(A) = 


BOF 


M(B, ) 


where B,, for i = 1,2, ..., k, are all the bases of [A:1I]. oO 
Before proceeding with our characterization of T(A), we should 

point out that a similar development is not possible for the integral 

property. The fundamental difficulty with the integral property 

is its connection with feasible bases. Which bases of [AjI] are 

feasible depends on the particular choice of b as a right hand side 

vector. Given an integer m yn matrix A, suppose we let Q(A) be 

the set of integer m-vectors b for which (3-2) and (3-3) define 

an integral convex polyhedron. Then, if pl and oe are in Q(A), 


bagi te F 2 
it is not true in general that bo +b is in Q(A). This is due 


70 


to the fact that the feasible bases associated with “ may not be 
the same as those associated with Bo Thus Q(A) is not closed 
under vector addition. 

The real advantage of dealing with the totally integral property, 
and our motivation in defining it, is that it involves all bases 
of [Ai I] for every choice of right hand side vector b. As a result 
it leads to an interesting algebraic characterization of T(A). 

We begin by considering a single basis B of [A:I] and showing 
how M(B) is characterized as the set of integer solutions to an 
equation expressed over a finite Abelian group whose structure 

L 


is completely known. An integer m-vector b is in M(B) iff Bb 


is an integer m-vector. If we let 
A = |det BI 
then BS is an integer m-vector iff 
AB +b = 0 mod A (3-5) 


where the columns of hace are integer m-vectors, 0 is the all zero 
m-vector, and the congruence applies componentwise. Consider the 
group MAB’), If we introduce the congruence relation on M(AB ¢) 
which relates two integer m=-vectors iff they are congruent modulo A 


componentwise, then the quotient group with respect to this congruence 
M(B”) 
MCAT) 


relation is » where I is the m xy m identity matrix. Moreover, 
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if we let Ba» Bo» sates BF denote respectively the congruence 
classes containing each of the columns of ABs and if we 
now let O denote the congruence class containing the all 


zero m=-vector, then (3-5) may be expressed as the following 


m(aB +) 


equation over the group G = M(AL) ” 


m 
_ b, B, = 0. 
tal 


Thus an integer m-vector b is in M(B) iff b satisfies the group 
equation (3-6) over G. 

We next show that (3-6) can be replaced by an equation over 
a new group which is isomorphic to G, whose size and structure 
is more easily revealed. First we note that if H is any 
group isomorphic to G, with isomorphism : G 4H, then the 
equation (3-6) over G is satisfied by some integer m-vector 


b iff the following equation over H 


m 
yb, o(B,) = 9(0) 


i=1 
is satisfied by b. 1 
M(AB ) 
We claim that the group G = -—“=—— is isomorphic to the 
M(aT) 


group G' The non-singular linear transformation - Bis 


= HM 
M(B) ° 


an isomorphism from M(B’) to M(I) and from M(aI) to M(B). It 


(3-6) 
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follows that ‘ B takes congruence classes in G into congruence 

classes in G' and that this mapping from G to G' is an isomorphism. 
We now take the isomorphism one step further. Smith [25] has 

shown that by performing certain elementary row and column operations 


on matrix B one obtains a diagonal matrix B, 


Ey 0 
Eo 
B = : 
0 a 


where €y> Foo see €, are positive integers, €,°€2° see Ge Ae 


and €5(s44 for i=1,2, ..., m-l. Matrix B may be expressed 


in terms of matrix B as follows, 


B=RBC 
where R-and C are integer m x m matrices with [det R| = [det C] = 1. 
We show that the group G' = neve is isomorphic to the group 
G" = aay - Since fdet R| = 1, R is an isomorphism from M(I) to M(I). 


Moreover, R is an isomorphism from M(B) to M(RB). But M(RB) = M(B), 
which can be shown as follows. Let y € M(B). Then y = RBCx = RB(Cx), 
for some integer m-vector x. Thus y ¢€ M(RB) and M(B) ¢ M(RB). Next 
let y ¢€ M(RB). Then y = RBx = RBC(C” ox) = B(C x), for some integer 
m=vector x. Since {det Cl = 1, oly is an integer m-vector and 


therefore y ¢ M(B). Thus M(RB) c M(B), and M(RB) = M(B). It follows 
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that R is an isomorphism from M(B) to M(B). Furthermore it follows 
that R takes congruence classes in G' into congruence classes in G", 


and this mapping from G' to G" is an isomorphism. 


M(I) 
M(B) 


the group of integers {0, 1, 2, ..., k-1} under addition modulo 


The structure of G" = is apparent. If we let J, denote 


k 
k, then it is readily seen that G" is isomorphic to a direct product 


group H, where 


The isomorphism from G" to H maps a congruence class in G" into the 
integer m-tuple obtained by taking any representative member of 
the class and reducing its {* component modulo Ey» for i=1, 2, ..., m. 
The order of H is Sy ee Although we have expressed 
H as an m-fold direct product, in many cases several of the E,'s have 
value 1, in which case the corresponding components of the direct 
product contribute nothing to the structure of H and can be 
omitted. 

The following diagram summarizes the sequence of isomorphic 


groups we have been discussing. 


1 reduction 
A B R modulo (E> Ego ses ED 
af —» —— —_—» 
NOE 2 cast, MO las, HCY: oe 
MTA: MGB “By. 7 ey, Peg: RR a 
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We began with the equation (3-6) over the group G as a means 
of charactérizing the members of M(B). We may now replace that 
equation with an equation over the group H. Letting o denote 


the isomorphism from G to H, we replace equation (3-6) with 


m 
A Nice, a0 (3-7) 
toe L L 
ic. 
where 
ey 
- € 
h. =o(B.) = 1 p(s 2 mod Bs 
A i . 
em 
aT 
€2 
= R, mod . 
L . 
Sai 


(The subscript i on a matrix denotes the — column of that matrix.) 
Equation (3-7) is the desired characterization of M(B). An integer 
m=vector b is in M(B) iff b satisfies equation (3-7). 

Now we are ready to characterize T(A), the group of all m-vectors 
b for which the constraints (3-2) and (3-3), with matrix A fixed, 
define a totally integral convex polyhedron. We recall from Lemma 3-3 
that 


k 
T(A) = TQ MCB.) 
i=l a 
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where B,, for i = 1, 2, ..., k, are all the bases of [Aj I]. 
For each basis Bs we have seen that M(B, ) may be characterized 


as the set of integer solutions to an equation 


over a finite Abelian group Hq, of order Ay = |det B.|- It follows 
that T(A) may be characterized as the set of integer solutions to 


k simultaneous equations, 
m 
> By hi, = 0, for i=l, 2, wees 
j=1 
over the finite Abelian groups H,, H,, sec, - Equivalently, 
1’ 2 


T(A) may be characterized as the set of integer solutions to a 


single equation, 


over a finite Abelian group K, where 


K =H, x Hy x tee x H 


and 


k. = (hy 5» Ros? eee, hy)» for j = 1, 2; eees Me 


The order of K.is A,‘A, © cee “Ale 
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We summarize these results in the following theorem. 


Theorem 3-1 


Given a fixed integer m y NM matrix A, the set T(A) of all 


integer m-vectors b for which the constraints 


[Afr] || =b 
Zz 


x 20, z20 


define a totally integral convex polyhedron is equal to the set of 


integer m=vectors b which satisfy an equation 


over a finite Abelian group K derived from all bases By > Bos eee, B.. 


of [A!I]. The order of K is A,*A5* ee "AL where 


A; = |det B.|- 


One conclusion we may draw from Theorem 3-1 is that T(A) never 
consists of just the all zero m-vector. If we let o(k;) denote the 


order of ks in the group K, and if we let @ be the diagonal m y m matrix, 
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0 


we see that 


M(@) < T(A). 


Thus T(A) contains a countably infinite number of mevectors. 

In terms of convex polyhedra having the totally integral property, 
these results are somewhat surprising. They tell us that for any 
matrix A there exist a countably infinite number of right hand 


side vectors b which give rise to totally integral polyhedra. 
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CHAPTER 4 


TRANSFORMATION OF (n,k)-POLYHEDRA INTO (k,k)-POLYHEDRA 


The general problem of integralizing convex (n,k)~polyhedra 
may be approached inductively with respect to one of the parameters 
nor k. Recall that an (n,k)-polyhedron is defined by Ax < b, where 
Ais anm yn matrix of rank kj, l<k <n. In this chapter we show 
that the class of (n,k)-polyhedra for fixed k and arbitrary n, 
n>k, have a great deal in common with respect to integralization. 
In particular we show that an (n,k)-polyhedron P,, where n 2k, may 


be transformed into a (k,k)=-polyhedron P, by a linear transformation 


k 


n k 
T: R »R which preserves integralization. 


In order to explain what we mean when we say that T preserves 


integralization, we refer to the diagram below. If we let ' denote 
T 
oo 
Po P. 
’ , 
p' ———__» P! 
n k 
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the operation of integralization, then this diagram expresses the 


fact that 
' = 1 
T(P)) = T(P) 7 


Thus in order to integralize Pa to obtain a we may instead integralize 
P. to obtain Pe and then use the inverse transformation ri to obtain 
Pe (A word of explanation concerning gt is in order. As we shall see, 
T is a many-to-one transformation. Thus it does not make sense to talk 
about tyy), ye R. as being a single real point x € Ry Rather, by 


they) we shall mean 
T (y) = (xeR | T(x) = y}. 


With this interpretation of cl we may say that tT (Pt) = Pre) 
Thus the purpose of this chapter is to show that the general 
problem of integralizing (n,k)-polyhedra reduces (somewhat) to the 
problem of integralizing (k,k)=-polyhedra. In the first section of 
this chapter we show how to obtain a suitable transformation T for 
a particular (n,k)-polyhedron Po In the second section we show 
that such a transformation preserves integralization. And in the 
third section we begin our inductive approach to integralization 


by examining the simple case, k = 1. 
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4.1 Existence of a Suitable Transformation 


Let Po be an (n,k)-polyhedron defined by Ax < b, where A is an 
integer m y n matrix of rank k, and b is an integer m=-vector. Any 
integer k yn matrix T of rank k whose rows span the row space of A, 
when viewed as a linear transformation from R” to aa transforms 
BP into a (k,k)-polyhedron P. defined by Cy < b, where y = Tx and 
C is anm y k matrix (not necessarily integer) of rank k defined 


by the matrix equation A = CT. However, we require a matrix T 


having two additional properties: 
(i) ged [T] =1 


(ii) every row of A is an integer combination 


or rows of T. 


In the following section we show that a matrix T having these 
properties preserves integralization. In this section we are 
concerned with the existence of such a matrix T. 

The conditions we have imposed on T are dependent only on 
the matrix A and are independent of the vector b. The following 


lemma assures us that a matrix T exists having the required properties. 
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Lemma 4-1 


If A is an integer m y n matrix of rank k then there exists an 
integer k y n matrix T of rank k such that (i) ged [T] =1 and 


(ii) every row of A is an integer combination of rows of T. 


Proof 


Let U be any k y n submatrix of A having rank k. According 


to Jacobson [18], U may be written as 


RUC 


Ga 
il 


where R is an integer k xy k matrix with det R = +1, C is an integer 


n yx nm matrix with det C =+ 1, and 
U = [u'; 0] 


where U' is a diagonal k y k matrix and 0 is the k y (n-k) zero matrix. 


Let T = [I ' 0], where I is the k y k identity matrix and let 
T=RTC. (4-1) 


It is clear that gcd [T] = 1. Jacobson [18] shows that gcd [T] = gcd [T] = 1. 
We now show that every row of A is an integer combination of rows 

of T. Let a be any row of A. Since A has rank k and Uis ak y n submatrix 

of A having rank k, the rows of U span the row space of A. The row a 


may be written uniquely as a linear combination of the rows of U, 
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i) 
W 

Q 

G 


where c is a rational k=vector (row). Then we may write 


a = cU = c(RUC) = cR(U'D)C = c(RU'R) (RIC) 


=c'T 


where c' is a rational k vector (row) given by c' = c(RU'R'!), 
We show that c' must be an integer k-vector. 

We may write c' = [c'J + <c'> where |c']| denotes 
the greatest integer function applied componentwise to c' and 
<c'> is the vector of fractional parts of c'. Each component 
of <c'> is a non-negative rational number less than l. 
Now, since a=c'T= [c'] T+ <c'>T and [c'] T is an 
integer n-vector, we see that <c's T is also an integer n-vector. 
We shall deduce from this that ec's = 0. 


We represent <c's uniquely by 
' i 
<C > = d (ey > C52 eee, e,) 


where d is a positive integer, C12 Cos coe, & are non-negative 


k 
integers, 0 < e,< d, and gcd (e,> Eas eres Ss d) =1. Let S 

be a subset of k colum indices of T such that the k y k submatrix 
T° comprised of the columns identified by S has non-zero 


determinant. Since <c'>T is an integer m=vector we see that 


Ss, , 
<c'>T is an integer k vector, say a'. Solving the system 


(4-2) 
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of equations a' = <c's r° for ec'> we find that <c'> has the form 


1 


TS are integers. 
|det T 


(£,> f5> eee fi)» where fj> fo. sais f 


Comparing this with (4-2) we conclude that d | det T°, Similarly, 


we may show that d divides all k x k minors of T. Thus d | gcd [T] = 


and d = 1. It follows that @, = ey Feee = Qe, = 0 and «c's = 0. 


4.2 The Transformation Preserves Integralization 


The transformation T whose existence is guaranteed by Lemma 
4-1 transforms the (n,k)-polyhedron Pi defined by Ax < b into the 


(k,k)-polyhedron Pe defined by Cy < b, where y = Tx and A = CT. 


Property (ii) of T insures that C is an integer matrix. We shall 


write P, = T(P) and Po - 17 (p,). In general, for a subset 


k 
k 1 . ‘ 
YoR, T  (Y) will have the meaning 


eye tx eR" | Tx ¢ Y}. 


We argue next that T takes faces of Pa into faces of Pe 


Let (Ps be an arbitrary face of Po 


Ps = {xe ea | 4g x = bg} 


1, 


O 
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and let (Pg be the face 
@MogtfyeR | Oy = bg}. 


Since A = CT we see that Ag - Co. It follows that TC(P)) = (Pigs 


It also follows that T*(@)g) = (Pp) The minimal faces of a 


s° 
are (n-k)-flats in R", Their images under T are minimal faces of 
Pe O-flats in R, 

Let ' denote the operation of integralization. Suppose we 
have computed Pao an integral (k,k)=polyhedron defined by C'y < b', 
where C' is an integer m' y k matrix of rank k and b' is an integer 
m'=vector. Then T (Pt) is an (n,k)-polyhedron defined by A'x < b', 
where A' is an integer m'y n matrix given by A' = C'T. The following 


lemma asserts that the transformation T preserves integralization. 


Lemma 4=2 


bw wielood 
Pt =T (Py) 


Proof 

We show that T (Pt) is integral and L(T"*(P!)) = I(P). It will 
follow from Lemma 1-1 that sage Se = Pa 

We first show that T* (Pt) is integral. Every minimal face of 
tet) is the image of a minimal face of P.' under cl and conversely. 


k 


Minimal faces of P' are vertices. Since P' 


k ” is integral, its vertices 
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0 
are integer points in rae Let y be any vertex of Pye Then 
ca Gra Oa ae al ee os Sa 
T Gy’) is a minimal face of T (Pi). T “(y ) is the (n-k)-flat 


defined by, 


thyy®s = {x @R" | Tx = ae 


Here we make use of property (i) of T, namely ecd [T] = 1. 


It follows that 
. 0 
ecd [T] = ged [Tiy ] 


and, by Theorem 2-1, vty?) is an integral flat. It follows that 
every minimal face of Ts) is integral, and by Corollary 2-4, 
Ty) is integral. 

We now show that L(t (Pty) = I(P). Since Py c Pe and 
coe) = P| we have mee) <P and L(t *(B!)) < IP). 
Consider the set difference Ee) - L(t" (Pt)). Assume there 
is an integer point x in this set difference. Since x ¢ I(P)> 
Tx = y is an integer point in P|, hence in Pe But then 
XE 1(T"*(P!)) which contradicts our assumption. We conclude 
that 1(P,) - L(t (P!)) is empty. Thus L(T*(P4)) =F: 

Having shown that T*(Pt) is integral and that L(t" (P!)) = 
1(P,), it follows from Lemma 1-1 that T"'(P!) = P'. 4 


Summarizing these results, we have seen that given a system 


of linear inequalities Ax < b defining an (n,k)=-polyhedron Pe 
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there exists a transformation T which transforms Pa into the 


(k,k)-polyhedron P, defined by Cy < b, where y = Tx and A = CT. 


k 
Furthermore, if we integralize P. to obtain Py defined by 
C'y < b', then Mg defined by A'x < b', where A' = C'T, 
is equal to pte Thus, without loss of generality, we may 
restrict ourselves to the problem of integralizing 
(k,k)-polyhedra. 

We examine the case k=l in the following section to 


illustrate the transformation and integralization of 


(n-1)-polyhedra. 


Awd Transformation and Integralization of (n,1)Polyhedra 


Let Pe be an (n,l)-polyhedron defined by Ax < b, where A 


is an integer m y n matrix of rank 1, and b is an integer m-vector. 


Let T = ass where a, is any non-zero row of A and A, = gcd [aj J. 


2s 
a 
Then gcd [T] = 1. Moreover, every row ae GE Ay doe Ly 2. ee yan, 
is an integer multiple of T, a, =, T, where Cc, is an integer equal 
to + ged [a, ]. Thus T satisfies the conditions for an integralization 


preserving transformation. 


IE) is the (1,1)-polyhedron P, defined by Cy < b, where 


1 
y = Tx and C is the integer m y 1 matrix defined by A = CT. 


According to Lemma 4-2, T*(P1) = Pie 
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The problem of integralizing Py is easily solved. Consider 
a single half-line 7; defined by c.y < by» Cc, #0. We integralize 


T, a8 follows. Let 1; be the half-line defined by 


b 
me F 
v= | ct ifc, > 0 


L 


Pe 
“y< Le | if Cc, < 0. 


The boundary point of Ty is an integer so that, by Theorem 2-2, 


Ts is integral. Furthermore, I(t,) = I(r,)- It follows from 
- is = t= 
Lemma 1-1 that Ti Ty H(I(t, )). 


Among all those inequalities cy < b, for which c, > 0, if 
b 
= is 
i 
smallest. Similarly, among all those inequalities c.y< b, for 


any, let the hi inequality be one for which the ratio 


which Cc, < 0, if any, let the ms inequality be one for which 


b, 
the ratio oa is smallest. Let Pi be defined by the following 
i 


system of at least one and at most two inequalities, 
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We see that 1) = I(P,). If - | > | then 
J 
1(P,) = I(P,) = ©, where @ denotes the empty set. Thus 
By = Pi =. Otherwise, by Corollary 2-4, By is integral. 
It follows from Lemma 1-1 that Py = Pie 
1 *(P!) is an (n,1)-polyhedron defined by the following 


system of at least one and at most two inequalities, 


o ca NOE ee 
By Lemma 4-2, T (Pi) = a 


Example 4-1 Let Py be the (3,1)-polyhedron defined by 


-8 4 12 xy 13 
12 -6 -18 Xo < 7 
“4 2 6 X3 9 
Let T = - [-8 4 12] = [#2 1 3]. Then Py = T(P,) is the 
(1,1)-polyhedron defined by, 
4 13 


89 


We integralize P, and obtain the (1,1)-polyhedron Ps defined by, 


1 


ys ‘ 
-1 1 


“2 1 3 3 
4 < 
2 = 3 1 
=3 


aL 1 = ' 
and T (P,) - Poe 


for i=1, 2. Since I(r;) = I(T.) we have 
= = ‘ ' 
1(P) = I(T, 9 T,) I(r) nN T)) 


so that P may be replaced by TT al Toe 


In the following we assume that this first integralization step 
has been performed on P. That is, we assume that P is defined by 


Ax < b where 


and 
gcd(a,y; 219) =1 
gcd(a,,s ayo) =1. 


Let A = det A. We assume that the order of the two inequalities 
Ax < b is such that A > 0. 
Our goal is to generate the system of inequalities A'x < b' 
defining P'. We accomplish this by generating all normal vectors 
and vertices of I(P). In the following we define the notion of 
a normal vector of I(P) and a vertex of I(P). We then introduce 
a transformation of co-ordinates which makes the properties of 
normal vectors more accessible. In particular, it allows us to identity 


images of normal vectors as being atoms of a certain partially ordered set. 
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We then develop a two-dimensional generalization of the division 
theorem for integers and show how it is used to generate the 
finite set of these atoms. We then describe a second level 
two-dimensional division process involving these atoms which, 
together with the inverse transformation, generates all normal 
vectors and vertices of I(P). These lead directly to the system 


A'x < b' defining P'. 


5.1 Normal Vectors and Vertices 
0 
0 as 
Let x = x0 be an integer point in I(P). We shall say 
2 


that - is a boundary point of I(P) iff there exists a non-zero 


oO 
integer vector g = ( *) such that 
02 


Se ay < 0 for all x € I(P) 


where . denotes vector inner product. We shall say that o 
is a support vector of I(P) at yo. The inequality ovteox’) <0 
defines a supporting half-plane for I(P). 

We shall say that a boundary point of I(P) is a vertex of I(P) 
iff there exist two linearly independent support vectors of I(P) 


0 
at x. 
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It is a straightforward matter to compute some of the boundary 
points of I(P). To do so we make use of some basic results in number 
theory. 


a 
Let m be an integral line defined by asx = b, where a = (2) 
2 


is an integer vector, b is an integer, and gcd (a> ay) = 1. 
An integer solution =" to a*x = b may be determined as follows. 


The Euclidean algorithm is used to compute a (non-unique) integer 


r 
0 
vector r = i) such that aer = 1. Then x = br is an integer 
2 


solution to a*x = b, and a-x = b may be written as eee) = 0. 
The following lemma is stated without proof. (For a proof 
see, for example, Niven and Zuckerman [21].) For convenience 


in stating it, we let 


Lemma 5-1 
If - is a line defined by han) = 0 where gcd (a,> a,) = 1, 
then I(m) is the set of integer points expressible as x + kTa, 


where k is an integer. 


0 


Consider the integral boundary lines 1, and T, of the integral 


1 2 


half-planes T, and T, defining P. Ty and TT, are defined by al.x =b 


1 2 1 
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a ag 


2 F = 
and a *x = by respectively, where a = ce and ged (a.7> a,») = 1, 


1 
for i = 1,2. The Euclidean algorithm gives two integer points x and 
2 
x in Wy and 7 


be written as a+ (x-x") = 0, fori =1,2. Lemma 5-1 tells us that 


> respectively. The defining equation for 7, may then 
I(m,) is the set of integer points expressible as so k,Ta", where 
ke is an integer. 
eae ; ‘ F 2 
Certain integer points in I(n) satisfy a-x< by and thus are 


integer points in I(P). It is easy to show that 1(7,) 1. ICP) is 


is 


pa where ky 


: F F 1 
the set of integer points expressible as x -k 


b Sa 
— i eee 
me ae 


Similarly, T(,) . I(P) is the set of integer points expressible 


1 


an integer and 


2 2 ; A 
as x + k,Ta , where ky is an integer and 


2 
bona bexe 
k fe eee 


Both TGr,) . I(P) and I(n,) 9. I(P) are sets of boundary points of I(P). 


0 
If x is a boundary point of I(P) then we shall say that a non-zero 


; é 0. 
integer vector g is a normal vector of I(P) at x iff: 
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(i) BGES”) < 0 for all x ¢ I(P) ) 
| 
(iis ar RTS: ae (PD ? (5-1) 
| 
LiL) ped (ay » Oy) = 1. } 


0 2 - 
The inequality %(x-x) < 0 defines a half-plane 7+ which is a supporting 
0 
half-plane for I(P). The boundary line 7 defined by o.(x-x ) = 0 


is a supporting line for I(P). 7 


contains at least two integer 

: 0 : Oe 2 sis 0 : ; Pao 
points, x and cither x + To or x - To, in I(P). The condition 
gcd (oy > Sy) = 1 is included to provide a canonical form for oa. 
If this condition were not imposed, then every positive integer 
scalar multiple ko of a normal vector « would also be a normal 
vector, and as we shall see, this would be inconvenient for our 
purposes. 

Two normal vectors of I(P) are evident. For every boundary 

es or ' Te 1 

point x in Ir) - I(P), a is a normal vector of I(P) atx. 


2 
Similarly, for every boundary point ee in I(n) 1 I(P), a is 


a normal vector of I(P) at gee 

Since normal vectors are support vectors, we see that a boundary 
point ae of I(P) is a vertex of I(P) if there exist two linearly 
independent normal vectors of I(P) at ae 

In the remainder of this chapter we develop a method for 


generating the system of inequalities defining H(I(P)). The method 
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involves generating all normal vectors and vertices of I(P). 
Clearly, if mY is a vertex of I(P) and og is a normal vector 
of I(P) at ee then the inequality actee < O defines a 
half-plane whose boundary line contains an edge of H(I(P)). 
It follows that the system of all inequalities of the form 
gintexe) < 0, where - is a vertex of I(P) and o is a normal 
vector of I(P) at er defines H(I(P)). 

We begin by introducing a transformation of co-ordinates 


which makes the properties of normal vectors more accessible. 


5.2 A_Hull-Formation Preserving Transformation 


Let us regard the matrix A of the system Ax < b defining 

2 2 

P as a linear transformation A: R2 >» R. For allxeR, 
A: X4y where y = Ax. Since A #0, A is a nonsingular 


; ‘ F 4 -1 
linear transformation. Thus the inverse transformation A 


exists. es is given by 


A is not an integer matrix in general. 
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We shall denote the image A(X) of a subset Xc R under A by e 


so that 
E F 2 
X =fyecR | y = Ax, x ¢€ X}. 


Since 


, 2 
Ps fee kh Ae bf 
it follows that 
2 an of D 
P =.py eR faye By 


Letting M(A) denote the set of all integer tembinations of columns 


of A we see that 
I(P) = {y € M(A) | y < b}. 


We claim that the transformation A preserves the operation 
of forming the convex hull of I(P). This is illustrated in the 


following diagram. 


| | 


A * 
HE CP) ) ae BCP) ) 
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This diagram expresses the fact that 
H(Z(P)) = H((P)") 


which we prove in the following lemma. 


Lemma 5-2 


HCL(B))” = H(T(P)*) 


Proof 
We first show that H(I(P)) c H(I(P) ). Let y € H(L(P)) . 


Then 


y = ACA < +r a +X ae) 
ai 0 1 2 


where \., 41, 4, 2 0, No +r + No = 1, and a ere a &-T(P).. 


es Same) 
We then have 


0 1 2 


2 he + ny + hoy 


) Le 


¢) 1 2 peers : “ 
where y , y, ¥ € [C(P) . Thus y € H(ICP) ). 


We next show that H(I(P) ) ¢ H(I(P)) . Let y € H(I(P) ). 


Then 
se 1 2 
yo May + MY + Noy 
oe a : O- 2% * 
where Noe Mao No 2 0, Ng ty + ro =1l, andy ,y,y e1(P). 
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We then have 


_ 0 1 2 
y= Ao Ax + Aj Ax + A, Ax 


1 2 
x + Nox ) 


0 
AQ ox + ry 


where ao x - € I(P). Thus y € H(I(P))*. 


We should point out that unless A = +1, A does not preserve 
the operation of integralization. That is, (p'y* # (P")', where 
' denotes the operation of integralization. In order to see this, 
consider P*, Pp" is a (2,2)-corner polyhedron defined by y < b. 
Since the vertex y = b of P is an integer point, p* is integral. 
That is, P* = (PY) ', However, if A & +1, then it is not true 
in general that the vertex x = ay of P is an integer point, 


in which case P is not integral. That is P' # P. It follows 


* * * 
that (P') #P =(P )'. 
Example 5-1 Let P be the (2,2)*corner polyhedron defined by 


23 


-3 1 “5 


Fig. 5-1(a) represents portions of P, I(P) and H(I(P)). Fig. 5~1(b) 


* * x 
represents portions of P , I(P) and H(I(P) ). 


O 
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It can be shown that the linear transformation A preserves 
the face structure of H(I(P)). The image of a face of H(I(P)) 
is a face of H(I(P)") and conversely. 

Let a be a boundary point of I(P) and let o be a normal 
vector of I(P) at ao Then egy y < 0 for all x € I(?). In 


* 
terms of I(P) we have 


“1 0 
oA “(yry ) < 0 G-2) 
* 0 0 * 
for all y € I(P) , where y = Ax is a boundary point of I(P) . 


Multiplying (5-2) by As 0 to clear fractions we obtain 
o-(AA™*) (yey) < 0 

or in terms of matrix multiplication only, 
o (AA™*) (y-y®) < 0. 


Letting 


y =o (al) (5-3) 


we have 


0 
y-Gyvr-y ) < 0 


* 
for all y € I(P) . We shall refer to y as a normal vector of L(P). 
0 
aty. 
* 
Since normal vectors of I(P) play a key role in our development 


we shall examine their properties in some detail. First of all, 
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we see that a normal vector y of T(P)” is related to the corresponding 
normal vector o of I(P) by equation (5-3), or equivalently, taking 


the transpose of each side, by 
-1,T 
YS AK os 


Letting Ae Gay! we have 


Letting M(A) denote the group of all integer combinations of columns 
of A under vector addition, we see from equation (5-4) that y € M(A). 
The defining conditions (i), (ii), and (iii) in (5-1) for a 

normal vector « of I(P) at ne may be recast as defining conditions 
for a normal vector y of L(P) at a as follows. If ig” is a boundary 
point of 1(P)" then a non-zero vector y € M(A) is a normal vector 


* 0) 
of I(P) at y iff: 
; a) * 
(i) ve (yry ) < 0 for all y E I(P) 


Gi) yo + tye 10)" 


(iii) Ny¢M(A) for all real A, Ocrhel. 


G4) 


(55) 


103 


We have already seen that condition (5=5)(i) follows from 


(51) (i). Condition (5-1) (ii) tells us that 


0. ato € 1(P)”. (5-6) 


Now AT = TA, as is easy to verify. Thus ATo = TAo = Ty and (5-6) 
becomes condition (5-5)(ii). Finally (5-1) (iii) implies that A° 
is not an integer vector for all realA, O<A< 1. Condition 
(5-5) (iii) follows. Thus, if go is an integer vector which 
satisfies (5-1), then y = Ao is an integer vector which 


satisfies (5-5). The converse argument can also be made. 


Example 5-2 Let P be the (2,2)-convex polyhedron defined in 
Example 5-1. Vertices and associated normal vectors of I(P) 


* 
and I(P) are tabulated below. 


vertices 


normal 


vectors 


vertices 


normal 


vectors 
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* 
Notice that tabulated normal vectors of I(P) are images of the 


1 
corresponding tabulated normal vectors of I(P) under A = 
3 


5.3 Normal Vectors as Atoms 


* 
In this section we show that normal vectors of I(P) are 
necessarily atoms of a certain partially ordered subset of 


M(A). We begin with the following lemma. 


Lemma 5-3 


* 
If y is a normal vector of I(P) at - then y 2 0. 
Proof 
According to (55)(i), 


0 * 
ye(yry ) <0 for all ye I(P). 


Recall that r(P)” is defined by 
* 
I(P) = fy € M(A) | y <b}. 


* 
Since 7 € I(P) , y e M(A) and y < b. Consider yt = y° = GE 


Since aia) is an integer vector, (4) ¢ M(A). Thus yt e€ M(A). 


(5-7) 
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* 
Also, since oo <b andaA>Qd, ye <b. Therefore yt € ICP). 


Setting y = ve in (5-7) we obtain “Yy A <90or 
¥, 2 0. 
ee ‘ 1 0 0 
A similar argument using y =y - A reveals that 


Y) 20. 


A normal vector y is, by definition, a member of M(A). 


If we let M(A)” be defined by 
ae ~ 
M(A) = {z € M(A) | z 2 0} 


then Lemma 5-3 tells us that y is a member of M(A)*. As we see next, 
a much stronger statement about y can be made. 

Consider the partially ordered set mat, <)- The ordering 
relation < is, as usual, taken componentwise. The element 0 € M(A)* 
is the least element or universal lower bound. 

An atom of (M(A), <) is an element a € M(A) such that a # 0 


and 

= + 
for all z € M(A). 
Theorem 5-1 


* 
If y is a normal vector of I(P) at x then y is an atom of 


(m(A)*, <). 
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Proof 
Let y be a normal vector of I(P)~ at We By definition, 
- ~ + 

y € M(A) and y #0. By Lemma 5-3 y € M(A) . 


Assume there is a z € M(A)* distinct from O and y such 


: 
that z < y. According to (5=5)(ii), 3° +TYE (P)*, We show 

that i + Tz and io + T(y-z) are also members of r(p)” and use 

(5-5) (i) to conclude that y.(4Tz) = 0. It will follow that 

z=Ay, for some real X\, O<¢ X < 1, which contradicts (5-5) (iii). 

We first show that 0 cae ae a we Since z € M(A) 
and TA = AT, it follows that Tz € M(A). Since re € M(A), 
oe + Tz € M(A). Now since yl +tTve 1(P)” we have a +Ty <b. 
Then since 0 << z< y we have - + Tz <b. Thus re +7z2 € r(p)*, 

We next show that jo + T(y-z) ¢€ L(P) Since (y-z) ¢€ M(A) 
and TA = AT, it follows that T(y-z) € M(A). Since y € M(A) we 
have y° + T(y-z) ¢ M(A). Again a +Ty <b. Since 0 < (y-z) < y 
we have © Pade) <b. Thus oe ee ert = 1(p)", 

Now we argue that y-(4Tz) = 0. Setting y = Me + Tz in (5-5) (i) 
we obtain y-GTz) < 0. Setting y = ie + T(y7z) in (5-5) (i) we 
obtain ye (4T(y~-z)) < 0. Since y-(+Ty) = 0, this gives -y-(4Tz) < 0. 
Thus y*(4Tz) = 0. 

Finally y-(Tz) = 0 implies that z = Ay for some real number i. 
Furthermore 0 <¢z<y, z2#0andz#y, imply that O<Ne< 1. But 


z=ay, OeNe 1, contradicts (5-5) (iii), 


O 
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Theorem 5-1 tells us that normal vectors of 1(P)” are necessarily 
atoms of (mca), <).- As we shall see, there are at most A + 1 such 
atoms. In the remaining sections we show how to generate these 
atoms and then show how they are used to compute H(I(P)*), and 


thus H(I(P)). 


5.4 A Two-Dimensional Division Theorem 


* 
Having shown that normal vectors of I(P) are atoms of (ma), <), 
our next objective is to generate these atoms. To this end we present 
the following theorem, which may be viewed as a generalization of 


the division theorem for integers. It applies to integer 2=vectors. 


1 


a 
The notation ( 
ag 


> b 
) < je) cael ze by and a, < by. All quantities 
2 


are integers. 


Theorem 5=~2 


al say 
Given ( ) and ( ) such that 
NE V9 


(2) > (‘2) > 0 
O< Vy < Vo 


u 
there exist unique q and (,°) such that 
3 
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and 
Oe ae 
O< Wo/ < Wy. 
Proof 
"3 
We first show that q and 6 \ exist. 
a/ 
Let q and Uy be determined by the division theorem applied to uy 
and Uys 
Uy = Uy 7 Uy, U, > Us 2 0. 
Since Uy > Uy > O, it follows that q > 1. 
Then V3 is given by 
V3 = Vy 7 Vy 
and since q>1 andO«< V1 < Vo it follows that 0< Vo < V3° 
"3 
We next show that q and € ) are unique. 
3 
Assume 
e eae i) 
ANG Wo / V3 


(5-8) 
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and 


where 


Uy. > fae 
Oe ee — i) ae 


Subtracting (5-9) from (5-8) gives 


1 
u u u 
2 3 a: 
ee ne 
Vo! V3 \v3 
Now from (5-10), 
Uy > u, > 0 
and 
O02 “U4 > TUy + 
Adding (5#12) and (5-13) gives 


- ' - 
uy > ur, u, > TU, 


which, together with the following equation from (5-11) , 


eat = eau 
(Q-q) uw = usr uy 


implies 


(3-9) 


(5~10) 


(5-11) 


(5-12) 


(5-13) 


110 


This, together with (5-11) implies 


O)- C9, 


DO 
“Y 2 
Given two ordered pairs of integers 2 ) and & ) such that 
1 2 
ba | > Uy >» 0 
Oz ee < . successive application of Theorem 5=2 yields 
a sequence of ordered pairs 
() > (2) > o) Pe ite ‘ee a Gs: @ @ = 
O< \Wy/ < \y/ < \Wy/ < < Whoo! < Woy < 
(514) 


The sequence (5-14) terminates when oe 0, which occurs for some 
finite k since Up> Uns eee, Uy is a monotonically decreasing sequence 
and uy is a positive integer. 

There is some symmetry in the statement of Theorem 5-2 which 
allows us to apply the theorem in the reverse direction. If we 
display the ordering relations among the components of (?) and oO 
as follows, 


@ > 2) > 0 
O< Uz! < Wy 
we find that we have two ordered pairs to which Theorem 5-2 may be 


Vv 
applied producing unique q and (,.) such that 
1 
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(2) 4 (2-2) 


tee 
0 < \u,) 


u 
Example 5-3 Let i = ~) a w) = (3): Successive application 


of Theorem 5-2 gives the sequence 
('3) (8) (oa) (se) (en) (207) 
3 8/ \21/ \34/ \47/ \107/° 


Applying Theorem 5-2 in the reverse direction, we may extend this 


(70°) (1) C3) @) (ea) Ga) Gr) Gor) . 


The sequence of first components in (5-14) 


se ila eal eo, 


is one which would be produced by the Euclidean algorithm applied to 
wy and Uys The set of common divisors of uy and Uy is equal to 
the set of common divisors of u, and Us449 


In particular, if uy and uy are relatively prime, then so are uy and 


for i = 1, 2, 0, kel. 


Uio. for i = 1, 2; cee, k=l. 
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The sequence of second components in (5-14) 
Vig > Vynz > °° > Vz 20 


is one which would be produced by the Euclidean algorithm applied 


to VE and Vee? Similar remarks about the set of common divisors 


of Vv; and Vaud for i = 1, 2, ..., k-1 apply here. 


Any three consecutive ordered pairs in the sequence (5-14) 


are related by an equation, 


Fe 8G 


i Vee i42 


for i= 1, 2, ..., k-2. These equations may be used to express any 
ordered pair in the sequence (5-14) as an integer combination of 


any two consecutive ordered pairs in the sequence. 


Example 5-4 In the previous example, ‘ee may be expressed as an 


integer combination of on) and & as follows, 


16\ _ 3 ( 1 
Gt) 5 (47): 
Similarly, ‘Ge may be expressed as an integer combination of 


Co) and (7°) as follows, 


(sa) = oe wee (3): 
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5.5 Generation of Atoms 


In this section we show how to generate the set of atoms of 
cm(A)*, <). Before doing so, we establish some preliminary results 
concerning M(A). 


Recall that A is given by 


>I 
tl 


where ged(a,4> ay») 1 and gcd(ay,> Ayo) = 1. For simplicity 


in our notation, we let 


where 44 = Ayp2 My2 = ~Apy> Ay = “Ay9 ANd ayy = ayq+ Thus 


gcd(q44> O19) = 1 and ged (51° m9) = 1. 


Let ot = Gas and a = Gee - The Euclidean algorithm may 
12 22 


r r 
be used to compute an = ( ae) and an we = ( a such that ae = 1 
12 "22 


2 
and Sa = 1. They are not unique. Let 81 be an integer such that 
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1 2 
8, =a er mod A, Og &, < A 
and let Bo be an integer such that 
2 1 
S =a er mod A, O<« Bo < A 


where A = det A = det A>»O. 

Integers By and 8, are uniquely defined even though vectors 
= and - are not. To see that By is uniquely defined, let - 
and (x2)! be any two vector solutions x to the equation 7 +x = 1. 
We may write of ex = 1 as a Gee) = 0. From Lemma 5-1, 

(x2)! = - + To for some integer k. It follows that 
Puce = oe mod a. Thus By is uniquely defined by (5-15). 


A similar argument reveals that g, is uniquely defined by (5-16). 
2 


Lemma 5-4 


Proof 


2 2 
By 8 = @ or ) @ “ey mod A 


[Gt er") Ger) - Gaiety’) (ehete”) | mod A 


[ Le (e* te”) | mod A 


1 mod A 


(5-15) 


(5~16) 
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2 
Now, let (rt)! and (r )' be any two vector solutions to the 


2 
equations oe = 1 and q °x = 1 respectively. Let 


rt = (r})! + k,To” 
where 
2 
a -(r')'=g, 
_—L eee 
1 A 
and let 
- = (x2)! + k,Ta 
where 
1 2 
8470 o(r > 
ky = Fi : 


Then ot -r! =i, ae = 1, meee = 8) and wre = Bo° We define 


matrices Cc) and Cy as follows: 


Goo = Faq ane | “O19 
c= GC ™ : 
Go, "22 12 O44 
We observe that det Cy = 1 and det C, = 1. Let Gy, = Ac, and 
Gy = AC, Then 
A 81 1 0 
Gy, = and Gy = . 


(5-17) 


116 


Lemma 5-5 


M(A) = M(G,) = M(G,) 


Proof 
We prove that M(A) - M(G,). Let z € M(A). Then z = Ax for 


some integer vector x. It follows that z = ACC, c,)x Ac, (Cy *x), 


2 cl 
i” 1, and Cy x is an 


integer vector. Thus z ¢€ M(AC, ) = M(G,) and we have shown that 


where ee is an integer matrix because det C 
M(A) < M(G,). 


Now let z¢€ M(G,). Then z = G,x = (Ac, )x = A(C,x) for some 


1 


integer vector x. Since C\x is an integer vector, z € M(A). 


Thus M(G,) c M(A), and we have shown that M(A) = M(G,). 


A similar argument gives M(A) = M(Go)-+ 


Lemma 5-6 


ged (A, 84) =1 and ecd (A, 85) = Ls 


Proof 


We prove that gcd (A, 81) = 1. From (5-17) we see that 


1 1 
ee C M(G,). By Lemma 5-5, M(G,) = M(G,), so that ea € M(G,). 


Thus i) = Gx for some integer vector x. Th#s says that 
2 


L= x, 4+ X) 8) 


from which we conclude that gcd (A, 84)| 1. It follows that 


gcd (A, 81) = 1. 
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A similar proof shows that gcd (A, Bo) le 
CO) 
We are now ready to generate the set of atoms of the partially 
ordered set mca)”, <). Suppose we apply our two-dimensional 


division theorem successively to the columns of G That is, 


1° 


suppose we let 
26) OQ) 
Vy 0 V5 1, 
and apply Theorem 5=2 successively to obtain the finite sequence 


uu u u U, u = 0 
5 CD (2)2 (2) 2 2 CZ Ce)? (9 


Vino’ < Mind k 
where k => 2. We shall show that the set of integer vectors 
appearing in this sequence is precisely the set of atoms of 


(M(A)*, <<). 


For convenience we let W be the 2 y k matrix, 


Mo° Bs bs ny is oy ag 
3 ke2 Vw OME 
u, 


i ‘ . 
whose columns w' = 2), i=i1, 2, ..., k, are the integer vectors 


appearing in the sequence (5-18). Thus 


V3 Ven? Ven] 


ae ae "3 bs oy i Bgl ; 
ts I VE 


(5-18) 
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We first argue that every column - of W is a member of 
M(A)*. As we noted in the previous section, the equations 
used to derive the columns of W allow us to express every 
column we as an integer combination of the columns of G,° 
Thus ee € M(G,). From Lemma 5-5, wi ¢ M(A). Furthermore, 


since u, 2 0 and v, 2 0, w e M(A)*. 


We next show that 


1 0 


ae _ [ 8, rn ° 


Veal 


We begin by showing that Yay = 1. From Lemma 5-6, gcd (u,> uy) = 


Furthermore, as previously noted, gcd (u, 5 Uap = ged (uy, Uy)» 
for i=l, 2, ..., k=l. Thus ged (up u,) = 1. Since uF 0 
and Un > 0, it follows that Un = 1. 

We now show that MoS A by expressing zi as an integer 


combination of the columns of G,. The equations used to derive 


1 
W are as follows: 
1 2 
wo=qw o-w 
2 3 
w = dow -w 
k=3 k-2 k-1 
w qn 3 w 
k-2 k-1 k 
Ww 


= q..2” "~ We 


From the last of these equations we have 


(5-19) 


k k-1 k-2 
Wo= doo W -w . 


Since wu = O and Uwe > 1, we see that qn? = Uno? We then 
write (5-20) as 


kL kel _ k=2 
WS ae 2 ao 


Then using (5-21) and equations (5-19) we have 


( u u yw fe oo 
Ne 3 AH? 7 Mel “K-2 
7 ke2 k=3 
= De ay bas 
k wett i wh 
bi a i+1 
rae eee 
= UY vig 
>= aye = nw s 


From the last of these equations we see that 


= Ael- g,°0 . 


Thus Vv. = Ae 


(5-20) 


(5-21) 


(5-22) 
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Finally we argue that v 


= 7 x, A + g_x 
write a tnx G,x for some integer vector x. Thus is Bs ( 1 2 Ee 
2 
from which we conclude that 
Uny = 81% pn mod A. 
Multiplying both sides by By we obtain 
Boy = 889% 24 mod A. 
Using Lemma 5-4 and the fact that Uy = 1, we have 
Vipn1 = 8 mod A 
and this, together with 0 < Veal < A, implies Vint ~ 82° 
Thus we may write matrix W as 
1 43 U-2 | 
We= | ¢, i eee Gy |. 
%3 V_n2 


W may also be generated from right to left by starting with the 


columns of CG and successively applying Theorem 5-2 in the 


reverse direction. 


Example 5-5 Let 


103 ~46 


kel = 89° Since woe Ee M(G,) we may 


i: 
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Then A = det A = (103)(61) - (-27)(-46) = 5041. The Euclidean 


algorithm gives 


21(103) + 47(-46) = 1 


9(-27) + 4(61) = 1 


Then 
&, = 9(103) + 4(-46) mod 5041 
= 743 mod 5041 
and 
Bo = 21(-27) + 47(61) mod 5041 
= 2300 mod 5041. 
Thus 
5040 743 1 0 
G, = G, = 
1 0 1 5 2300 5041 


Matrix W is generated using either G, or G, and Theorem 5-2, 


5041 743 
2300 5041 


We require two additional lemmas before we prove that the set 


of columns of W is the set of atoms of (M(A)*, <)- 
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Lemma 5-7 
Let Gis i=1, 2, ..., k-1, be the 2 y 2 submatrix of W 
given by, 
"5 MEHL 
G. = . 
fs v 
i i+l 


Proof 

The equations used to generate the columns of W allow us to 
express the columns of G, as integer combinations of the columns 
of G,- Thus M(G, ) ¢ M(G,). The same equations allow us to 


express the columns of G, as integer combinations of the columns 


1 
of Ge Thus M(G,) ¢ M(G,), and we have shown that M(G,) = M(G,)- 


By Lemma 5-5, M(G,) = M(A). 


Furthermore, from equation (5=22) we see that 


MT Ga 7 aaa 


and since a A, we have det G, =A. 


We define Mt (Ww) to be the set of all non-negative integer 


combinations of columns of W, 


Mt (w) ={zeé qe | 2 = Wx, x € nae x > 0}. 
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Lemma 5-8 


Mt(w) = M(A)* 


Proof 

Since each column of W is a member of M(A)*, it follows that 
every non-negative integer combination of columns of W is a member 
of M(A)*. Thus M'(W) ¢ M(A)*. 

“+ + =\+ 

We now argue that M(A) « M (W). Let z € M(A). If z = 0 then 
clearly z € Mt (Ww). Therefore we can assume that z # 0, that is 
z4.> 0 or Zo > O or both. 


Consider the quantities 


for i=l, 2, ..., k. For i = 1 we have 


21 uy Z4 A 
= =-z, 4 <0 
0 
Zo vy Zo 
and for i = k we have 
Z4 uy. zy i¢) 
= =z, A 20. 


(5-23) 


(5-24) 
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We claim that 


1 i 1 0 Mad 


2 i o. “saa 


(5-25) 


for i=l, 2, ..., k-1. To prove this we observe that 


0 Vv. Vv. and u, u, 
Mi < ‘GG eat E 


+1 


> 0. 


(5-25) follows from this and our assumption that Z> 0 or Zy > 0 or both. 


Now, from (5-23), (5-24) and (5-25) we conclude that either: 


(i) for some i, Lleiek 


or 


(ii) for some i, leiek 


We consider case (i). By Lemma 5-7, M(G,) = M(A). Since z € M(A)* 


it follows that z ¢€ M(G,) and z = G,x for some integer vector x. We show 


that x, > O and x, > 0. From Lemma 5-7, det G = A. Applying Cramer's 
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rule, we solve z = G,x for x: 


24 Vert 
1 
=i 0 
xy A > 
49 Nig 
i oy a Uy 
xt Sag 
a aa “7 A ae 
Ve 29 a) Be 


Thus z is a non-negative integer combination of columns of G,- 


u, 
i 
xy e) where xy 


In case (ii), similar reasoning gives z 
is a positive integer. 
+ =\+ + 
Thus we have shown that z € M (W). It follows that M(A) cM (WW). 
+ =\+ 
We conclude that M (W) = M(A) . 


a) 


We are now ready to prove the main result of this section. 


Theorem 5=3 


= + 
The set of columns of W is the set of atoms of (M(A) , <). 


Proof 
‘ oak 
We first show that every column w of Wis an atom of (M(A) , <). 


‘ P aes 
We know that w- E M(A)* and w #0. Assume that there is a z € M(A) 
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distinct from 0 and w’ such that z <w. By Lemma 5-8, z may be 


expressed as 
k 
Zz »: c 7 
i 
i=l 


where Cy is a non-negative integer, for i =1, 2, ..., k. If c, = 0, 


for i=1, 2, ..., k, then z = 0 contrary to assumption. Thus 


c, > O for some j, 1 <¢ j < k, and it follows that w! < Zz. Since 


= 
tu. 


<zandze< we it follows that w < ws which implies that 


w = w', because 
| > Uo > eee au. 
Vy < YW <ttt< Ye 


Now we have w <zZzc¢ a and w as We conclude that z = w which 


again contradicts our assumption. Therefore no such z exists and 
w is an atom of (M(A)*, <)- 
We next show that every atom of (M(A)*, <) is a column of W. 


Let a be an atom of (M(A)*, <). Then a # O and 


Zea =>» 220 orze=a (5-27) 
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where cy is a non-negative integer, for i=l, 2, ..., k. Since 
a #0, there is some ey > 0, 1 < j < k, in this expression. 
It follows that w” <a. According to (5-27), w =Oorw =a. 


We know that w? #0. Thus a = w 


Matrix W has at most A +1 columns. This follows from the 


strict ordering, 


Theorem 5~3 tells us that (M(A)", <) has at most A + 1 atoms. 


5.6 Generation of Normal Vectors and Vertices 


We have seen that a normal vector y of 1(P)" at a vertex y 
is an atom of the partially ordered set (M(A)*, <), where A is 
the cofactor matrix of A and Ax < b defines P. We have developed 
a method for generating the matrix W whose columns are the atoms 
of (M(A)’, <)» In this section we complete our development of an 


integralization method for (2,2)-corner polyhedra by showing how 
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all normal vectors and vertices of 1(P)", hence of I(P), are generated. 
Our final product will be a system of inequalities A'x < b' which 
defines P' = H(I(P)). . 

We begin by considering the following question. Given a 
boundary point a of 1(P)", how do we tell which columns of W 
are normal vectors of r(p)” at yi? 


i * i * 
Let yo be a boundary point of I(P) . Since y" € I(P) we have 


a <b. Let z = bey. Then eS > 0. z° measures the "slack" 
of vs 
Recall that a normal vector y at yt satisfies the following 
conditions: 
i) ye (y-y") < 0 for all y € (P)* 


(ii) y' + Tye 102)". 


Using the slack vector zy we label two columns of W as shown 


below: 
A = uy = 0 
O0O=Vv =A 


- is the largest element in the lower row of W which is less than or 


equal to zy it is the largest element in the upper row of W which is 


less than or equal to 25% 
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Theorem 5-4 


P q 


i * 
Tf y" is a boundary point of I(P) then the columns w’ and w 


:< i z 
of W are normal vectors, Y" and yank respectively, of I(P)” at ue 


Proof 
* « 
According to Theorems 5-1 and 5-3, normal vectors of I(P) at y* 
are columns of W. 
‘ F itl i 
We first observe that there exists a normal vector y at y 
i i+l * ‘ 
such that y - Ty € I(P) . We shall consider all columns w of W 
i * A i+l 
such that y - Tw € I(P) to be candidates for y and show that 
i * 
all but w? fail to satisfy the condition we (yy) < 0 for all y € I(P). 
It will follow by elimination that yr =wi, 

For all columns w of W, y+ Tw € M(A). This follows from the 
following argument. Since w € M(A)*, w = Ax for some integer vector x. 
Then Tw = TAx = ATx, so that Tw ¢ M(A). Since y' ¢ M(A), y + Tw € M(A). 

i * i 
Thus y +TweI(P) iff y +Twsb. 


J 


Now, for all columns w’, where j <q, y- -twi £ b because 


Y> + u, > bos For all columns w, where qs j <k, y - Tw <b 


i i i j * 
because yy v5 < by and Yo + u, < bos Thus y' - Tw ¢€ I(P) 


for j =4q, qtl, ..., k. However, for all column w, j = qtl, qt2, ..., k, 


wi. (yt ~ Ww) -y4) = we (-twt) = $0 


oe : x é : —— a ais 
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because a > 4, >O0O and ws > ve 20. Thus all of these columns 


‘ * 
fail to satisfy the condition wie (yey?) < 0 for all y € I(P). 


Therefore, by elimination, yr = wi, 


iel i 
We next observe that there exists a normal vector y at y 


; ee x 
such that y- + Ty" ¢ I(P) . A proof similar to the one above 


establishes that or ge =w, 0 


In Theorem 5-4, the columns w’ and wi need not be distinct. 


; es ‘ 
If they are distinct, then y- is a vertex of I(P) . Otherwise y- 
* 
is merely a boundary point of I(P) . 
i * 
If y is a vertex of I(P) then Theorem 5-4 gives us two 


inequalities 


ane oe <0 


(5-28) 


oft}, (gexty <0 


in the system A'x < b' defining H(I(P)). In (5-28), Pi = ay" yett 
and xt = ae yi. 

We now consider the following question: given a vertex y* of L(P)*, 
how do we generate neighboring vertices, if any, of y'? 


. ‘x : 
Let y™ be a vertex of I(P) . We assume that z= beyt has been 


used to label two columns of W as before. We show how to generate 


slack vectors gee = bey? and a = bey of neighboring vertices 


i-2 i+2 F 
y and y respectively. 
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We refer to the previous labeling of W. If Le > 0 (that is, 


p > 1) then let 
. ‘ A 
go See Is, | tw, (529) 


Otherwise, disregard 2 If oS > 0 (that is, q < k) then let 


. ; Zal 
ae eee lr Tw, (5=30) 
q 


: L-2 2 
Otherwise, disregard oot Using z and gre we label two 


additional columns of W as shown below: 


i 
Oo 


v.. is the largest element in the lower row of W which is less than 


i~2 


or equal to 2; . ms is the largest element in the upper row of W 


which is less than or equal to “Y 


Lemma 5-9 


or ‘ a 
y- ‘ and yor? are vertices of I(P) . 
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Proof 


i+2 . * ; q 
We prove that y is a vertex of I(P) by showing that w 


* i42 
and w are two distinct normal vectors of I(P) at yr - We have 


, Zoi 
yi? “ yt = l= | tw". 
“q 


$42 P ei 
Clearly yo <b, so that ge c I(P) . Since wieTwi = 0, it follows 


; P i 
that et ea = wie (yey) < 0 for all y ¢ I(P) . Also, since 


z 
: ‘ * 
l= | 21, yor + Twt ¢ I(P) . Thus wi is a normal vector of 


* ‘ Pa 
I(P) at A aes By Theorem 5-4, w° is a normal vector of I(P) at 


i+2 . i+2 i42. 
y- - Since a0" < a we have we # wi. Therefore y- is a 
* 
vertex of I(P). 
ete i-2, * 
A similar argument reveals that y is a vertex of I(P). 
O 
What we have described may be viewed as a second level 
two-dimensional division process in which the slack vectors 
of known vertices are the dividends, atoms of (M(A)™, <) are 
divisors, and slack vectors of new vertices are the remainders. 
In order to initiate the process we must have an initial 
vertex. This is easily obtained. In Section 5.1 we saw that 
I(n) f ICP) and I(r) MN I(P) ave sets of boundary points of 


I(P), where mn, and T, are boundary lines of half-planes defining P. 


1 
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Furthermore we saw that T(r) f ICP) is the set of integer points 
expressible as x}ak tal, where = is an integer point in (7) 
computed using the Euclidean algorithm and ky is an integer such 


that, 


b ee at 
ky < || ; 


We obtained a similar characterization of I(t) fn ICP). 


tac -8 a aaa 1 
If we take x ~k,Ta » where ky = | » as being x , 


then we claim that is a vertex of I(P). Consider - = hee 


If we compute oe = Baye we find that z = 0 and 25 < Ae 


P and wi are normal vectors 


It follows from Theorem 5-4 that w 
* 1 1, 
of I(P) at y', where p =1 and q>1. Thus y is a vertex 
* 
of I(P). 
Similarly, we can find an initial vertex drawn from 
I(m,) 9 I(P). 
‘ 1 Pastas 1 * 
Using the slack vector z of the initial vertex y” of I(P) , 
we label two columns of the matrix W as we did for Theorem 5-4, 
The circled slack vector oe then points to the normal vectors 
0 2 * 
y andy of I(P) at _- The division indicated by equation 
(5-30) is then performed on the labeled matrix W to obtain slack 
vector - which in turn is used to label a third column of W. 
The circled slack vector oe then points to the normal vectors 


2 
yY and we — and so on. The diagram below illustrates the labeled 


matrix W. 
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Example 5-6 Suppose W is given by 


73 56 39 22 5 3 #1 =O 
0 1.2 3 4 17 30 73 


0 
70 


ve Then the labeled matrix W is shown below. 


5 3 1 0 
4 17 30 73 
Sample divisions are: 


(4) = (8) +2] 
() = (2) +P] (4)- 


Notice that Equation (5-29) enables us to check our division by 


ar 
and suppose z = 


working back from right to left. 


O 
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: 2n-1 . 
The process terminates when a vertex y is found such that 


2n-1 . eee 
Zo = 0, which must occur after at most k-1 divisions have been 


performed, where k is the number of columns of W. If the division 


process is executed from right to left, beginning with a slack 


-1 ; -2n+1 . 
vector z , then the process terminates when a vertex y is 


-2nt+1 = 


1 e 


found such that z 
In the final labeled matrix W, each circled slack vector 


* 
1 of IP) 


i F i-l i+ 
Zz points to the two normal vectors y and y 
i i. , 
at vertex y'. Each normal vector y is pointed to by the one 
* 
or two circled slack vectors representing vertices of I(P) 


at which Y- is a normal vector. 


The final result of all this is an alternating sequence 


* 
of normal vectors and vertices of I(P) , which may be transformed 


into the alternating sequence 


itl ,=,*1 i+l 
oo” Y a 


of normal vectors and vertices of I(P), where =(A) nd 


x = anty*, for i = 1, 3, ..., 2n-1l. We then write the following 


system of inequalities, 
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orien) < 0 
ees) <0 


tes) <0 


Since 2 = 0 and z = 0, it follows from Theorem 5-4 that 


v = @ and yen G)- Thus Ps = Gy"! ye = at and 


0 
oo = yoy = vt = b,, so that the first inequality 
in (5-31) is ioe < by: Similarly 28 = (A)7} ee = ae and 
gga’ = 4 yey ae = yon = b,, so that the last 


A Z 2 


inequality in (5-31) is ance < bo- These two inequalities 
are the two which define P. 

If we let Q denote the polyhedron defined by the inequalties 
(5-31), then it is easy to see that Q = H(I(P)). Since each of 
the half-planes whose intersection is Q contains I(P), it follows 
that I(P) ¢ I(Q). Since the two half-planes whose intersection is 
P are included, it follows that I1(Q) © I(P). Thus I(Q) = I(P). 
Furthermore, by construction, every vertex of Q is an integer point. 
Thus Q is integral, From Lemma 1-1, Q = H¢I(P)). 

The following lemma gives an upper bound on the number of 


vertices of H(I(P)). 


(5-31) 
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Lemma 5-10 


H(I(P)) has at most A vertices. 


Proof 
We have previously argued that matrix W has at most A + 1 
columns. Therefore, in labeling the columns of W, at most A 


vertices are generated. 


We conclude this chapter with two examples illustrating 


the integralization of (2,2)-corner polyhedra. 


Example 5-7 Let P be defined by 


[5 31 Ei <[7} 


Then 
is S 
a i ri 


138 


An initial vertex Pa of I(P) is given by, 
where 


0 


Thus x = Cae and a 


The resulting system of inequalities defining H(I(P)) is: 


(5) + (o)) = ° 


1) * oe = > : 


i) 


wo 


= (Gy) W is labeled as follows: 
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Example 5-8 Let P be defined by 


279 


s 


Xo -317 


Using the Euclidean algorithm we find that ged (488, 216) = 8 and 
gcd (46, 103) = 1. Since 84279 the first inequality does not 
define an integral half-plane. The second one does. We replace 


the first by 61 x1 + 27 X < |e and obtain a new polyhedron 


P, defined by 


1 
61 27 xy 34 
s 
46 103 Xo -317 
Then 
‘ 103 -46 
A= - 
=27 61 


Matrix W has been computed in Example 5-5. An initial vertex a 


of I(P,) is given by, 


-9 ~61 
where 
seeks ios -34( 4, 
ie wh 
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Thus a = Gy and so = Ge \ W is labeled as follows: 
vA 


977 1859 


1418 


The resulting system of inequalities defining H(I(P,) is: 
(7) + (Ce) = ° 
27 *"\462/)) § 
i) -(x- (62) < © 
\. (x- -(4 )) < 0 or 2 1 0 
17 38 -118 

3) + (e(4)) <0 

36 

03 (x -(33)) <9 


61 27 x 34 


tN 


46 103 ~317 


(: 
(r 
(3 
(ro 
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CHAPTER 6 


INTEGRALIZATION OF (2,2)-POLYHEDRA 


Let P be a (2,2)-polyhedron defined by Ax < b, where A is an 
integer m y 2 matrix of rank 2, and b is an integer m-vector. 
In this chapter we apply the integralization method of Chapter 5 
to the problem of integralizing P. 


We begin by reducing P to an intersection of integral half=planes 
F 14 by 
by replacing each inequality ax < b, in Ax < b by 4 ax < a, 
i i 
where d, = gcd (ajy> as o)s If we let Ts denote the half-plane 


defined by ax < b, and rj denote the half-plane defined by 


b 
1 i i ; P 
a, a-x< las >» then T(t;) = T(r;) for i=1, 2, ..., m. Since 
m m m 
1(P) =  1(t,) = M T(r{), P may be replaced by 1; without 
i=l i=l i=l 


changing I(P). 
In the following we assume that this first integralization step 
has already been performed on P. That is, we assume that each row 


£ = e 
las4 ao] of A has the property gcd (a,4> a.) = 1 


i2 
Our plan is to integralize P by generating all vertices and 


normal vectors of I(P). Boundary points, vertices and normal vectors 
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of I(P) are defined just as they were for (2,2)-corner polyhedra 
in the previous chapter. We accomplish this by generating 
boundary points and normal vectors for a sequence of (2,2)-corner 
polyhedra which we shall call supporting corners of I(P). 
A (2,2)-corner polyhedron Q is a supporting corner of I(P) 
iff I(P) ¢ I(Q) and I(Q) and I(P) have one or more vertices 
in common. 
In the following we show how an initial supporting corner 
of I(P) is found. We then show how supporting corners of I(P) 
are used to generate vertices, normal vectors and other supporting 
corners of I(P). We conclude by showing that the process terminates 
after a finite number of steps, resulting in a system of inequalities 


A'x < b' defining P'. 


6.1 Finding gn Initigl Supporting Corner 


Well-known methods associated with linear programming (see, for 
example, Dantzig [5]) may be used to find a vertex v of P. Associated 
with vertex v is a two-inequality subsystem Cx < d of the system 
Ax < b defining P, where v = cla, Cx < d defines a (2,2)-corner 


polyhedron Q such that I(P) ¢ I(Q). 
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Applying the methods of Chapter 5 to Q, we obtain the sequence 


x x eee 


of normal vectors and vertices of I(Q). Suppose x € I(P) for some i, 
; : i+] i 
where i = 1, 3, ..., 2n-l. Since o ™. (xx) < O for all x € I(Q), 
; itl i 
it follows from the fact that I(P) ¢ I1(Q) that go —.(x=x’) < O for 
all x € I(P). Thus x is a common vertex of both I(Q) and I(P), and Q 
is a supporting corner of I(P). 
Now suppose x" ¢ I(P) for i = 1, 3, ..., 2n-l. In this case, 
P either does or does not intersect the boundary of Q' = H(I(Q)). 
Suppose P does not intersect the boundary of Q'’. Then since 
the vertex v is in P and P is convex, it follows that P lies 
wholly within Q-Q' and thus I(P) =o and P' =, where o denotes 
the empty set. 
Suppose on the other hand that P does intersect the boundary 


of Q'. Let e, for i=l, 3, ..., 2n-3, denote the edge 


‘ P : : 
(x E Ro | x =x 4 iC coal tay Oe 21) 


of Q'. We argue that the intersection of P with the boundary 
of Q' is containing in exactly one edge a” of Q', as illustrated 


in Fig. 6-1. 


Fig. 6-1 


First, we see that neither the half-line edge 


{x € R? | x xt bite oy) X% = 0} 


nor the half-line edge 


{x € R? = ae + a ame -v), \ 2 0} 


of Q' contains any points in P. Otherwise, since v € P, the convexity 


2n-1 


of P implies that x € I(P) or x € I(P), contrary to assumption. 


Next we see that if two edges Qe, and e, of Q', where i < j , both 
contain points in P, then since v € P, the convexity of P implies that 
ae oo. ee x € I(P), again contrary to assumption. Thus the 
intersection of P with the boundary of Q' is contained in exactly 

' 


i 1 
one edge e” of Q'. We form a new (2,2)-polyhedron P by adding 


the inequality 6 eee) < 0 to the system Ax < b and deleting 
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the subsystem Cx < d. Since p-pt lies wholly within Q-Q' and thus 
contains no integer points, we see that (pt) = I(P). Thus we 
have reduced the number of inequalities defining P by one, 
without changing I(P). 

Computationally, these observations are implemented as follows. 
For each vertex xt of 1(Q), i = 1, 3, ..., 2n-1, we compute the 
slack vector ze = peda If oe 2 0 for some i, then xt € ICP) and, 
as we have seen, xt is a common vertex of both I(Q) and I(P). Thus 
Q is a supporting corner of I(P). 

Otherwise zt % O and xt ¢ I(P) for i =1, 3, ..., 2nel. 
In this case, for each i, i = 1, 3, ..., 2n-3, we form the following 


system of m + 2 inequalities in \, given below in vector form as, 
zo + nz teazty > 0 


O<AX <1. 


(6-1) has a solution \ iff P intersects the boundary of Q' at edge 
ae As we have seen, two cases are possible. In the first, (6-1) 
has no solution dX, for i = 1, 3, ..., 2n-3. In this case I(P) = 9. 
In the second case, (6-1) has a solution X\ for exactly one i among 
1, 3, ..., 2n=3. We then form a new (2,2)-polyhedron pl by adding 
the inequality ot, (xox?) <0 to Ax < b and deleting the subsystem 
Cx <d. 
If we have not found a supporting corner of I(P) or learned 


that I(P) =, we may apply the entire procedure described above 


(6-1) 
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to pt, In general, during the oo application of this procedure 
we either: (a) discover a supporting corner of I(P), or (b) find 
that I(P) =, or (c) obtain a (2,2)-polyhedron ps defined by 
m-k inequalities such that tie) = I(P). 

We terminate this procedure when either (a) or (b) results 
from a given pass. If neither (a) nor (b) have resulted from 
the first m-2 passes, then we are left with a (2,2)-corner 
polyhedron oo such that tPF?) = I(P). By definition, pmn2 
is a supporting corner of I(P). 

We assume now and for the remainder of this chapter that 
T(P) #. Thus a supporting corner of I(P) is found on the ge 
pass, where 1 <¢ k ¢ m-2. In the following, for convenience, we 


refer to the resulting (2,2)-polyhedron pol as P and to its 


defining system of inequalities as Ax < b. 


6.2 Generation of Normal Vectors and Vertices 


Let o* be a supporting corner of I(P) and let ra be a common 
vertex of both T(Q*) and I(P). We refer to the pair (aK, 5 as a 
State of the integralization process. gf is defined by two 
inequalities ck < ats We assume they are ordered such that 


det ck > 0. We shall refer to the inequality toward whose boundary 


line we are moving as the objective inequality. The objective 
F : : : k , 
inequality associated with Q° is drawn from the system Ax < b 
defining P. We shall denote the objective inequality associated 
; k J th. 
with Q by a «x < » the jy inequality in Ax <b. 
k 


We proceed from state to state by means of a basic iterative 
step. A state which results from one application of this iterative 
k i 5 F : k k 
step to (Q’, x), where the second (first) inequality of Cx «< d 
. . : . , . k i 
is the objective inequality, is A+ (-) successor state of (Q°, x). 


kk : 
The + successor of (QQ, x") takes one of the three forms illustrated 


in the diagram below. 


Cis) 
eee | Ps 
ee kt] _i+2, eel st 


(Q*, _ ) (Q air a? ae 2 (Q 


k ne : 
The + successor of (Q , x ) is characterized by a new vertex, 


a new supporting corner, or both. 


Example 6-1 Fig. 6-2(a) and (b) illustrate a succession of states: 
k i : , , 
COr A ear SOO ee” Se SE ia a a CO ae hg 
C 
0 0 
A state (Q , x ) used to start the process is an initial state, 


0, , . ) 
Q is defined by two inequalities C x < a both drawn from Ax < b. 
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Thus either inequality may be designated as the objective inequality 


Jo 
a ex <b, . We saw how to find an initial state in the previous 
Jo 
section. 
mF 
OBJECTIVE OBIEcTiwve soe 
INBQUALITY INEQUALITY 
(a) (b) 
Fig. 6-2 


A state which has no + successor or no - successor is a terminal 


state. 


We now describe the basic iterative step for determining the 


+ successor! of (Q, x’), if one exists. 


t For + successor, all instances of +(+) should be read as +(-). The 


opposite interpretation holds for = successor. 
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(1) Using the (2,2)-corner polyhedron integralization method 


of Chapter 5, we compute the normal vector gith of 1(Q*) at x. 


(2) We let 6x? = Trott! We compute 


Sy> = A 6x" 


and compare by" with the slack vector zy where 
z = beAx” = bey. 


(3) If all components of 6y- are non~positive then x +n 6x € I(P) 
: itl. i 
for all integers n 20. o — is a normal vector of I(P) at x”. We 
adjoin the inequality Pe Om < 0 to an evolving system of inequalities 
defining P'. The state (QX, x) is a terminal state and has no + 


successor. 


(4) If one or more components of $y" are positive then we compute 


the non-negative rational numbers 


for all subscripts j such that by, > 0, and among these we identify 


a subscript Jad such that 


Nici <5 
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for all j such that by,” > 0.7 We distinquish three possible cases 


corresponding to the three possible forms of the + successor of 
k Ui é og , . ; : A 
(Q°, #0 (a) jg = Jy (b) Jeg # j,, and rey 21, (c) Jeg Fy 


and 0 < Niad <i. 


: er : F 5 i i 
(5a) Get - j,) In this case NI, 21. Since x” + 6x” € I(P), 

ot is a normal vector of I(P) at x". We adjoin the inequality 

7 ; 
one (x=x") < 0 to an evolving system of inequalities defining P'. 
We let 

oe =x + |x| 8x". 
k 
i+2 


x“ is a vertex of both r(Q°) and I(P). We take ck, yrs as the 


+ successor of cok, x). 


(5b) Gres # ii and Mice 21) In this case, since 
x + 8x e i@); Papa is a normal vector of I(P) at x. We adjoin 
the inequality git tae) < 0 to an evolving system of inequalities 


defining P'. We let 


ae & . a 
x =x + | Sect | x. 


The following somewhat arbitrary rule insures the determinism of 

this step: If there is more than one subscript j such that Aj 

is smallest, then choose i, if it is included — otherwise choose 
the smallest subscript. (Recall that i, identifies the objective 


inequality associated with g*.) 
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k+l 


In addition, we let Q be the (2,2)-corner polyhedron defined by 


; ; j 
ott. (ext) < 0 and the new objective inequality a Kil < bia 


i+2 
drawn from Ax < b. We see that I(P) ¢ rt!) and that ait 


is a vertex of both rq!) and I(P). Thus htt is a supporting 


corner of I(P). We take (ght | xi?) as the + successor of 


ca, x4). 


(5c) Grey # j,, and O< Nia <1) In this case, since 
i i : F k+1 
x + 6x” ¢ ICP), no new vertex is obtained. We let Q-— be 
the (2,2)-corner polyhedron defined by tec’) < O and 
j 
the new objective inequality a Kil < bid drawn from Ax < b. 
We see that I(P) ¢ rcg'ty and that x’ is a vertex of both 


ratty and I(P). We take ght, x") as the + successor of 


k Ji 
(Q’, x). 
This completes our description of the iterative step. 


The integralization process proceeds as follows. Starting 
P an eg he 0 0 + 
with an initial state (Q., x ), a sequence S of + successors 


- F k 2n, , . 
is generated until a state (Q., x ) is obtained such that: 


(a) n>2 and xen - = 


or 


(b) cK, ny is a terminal state. 


ki 2m . , ; : Peete 
If (Q°, x) is a terminal state then starting with the initial state 


0 _0 - F : 
(Q., x ), a Sequence S of = successors is generated until a terminal 


wl! =2n'! 
state (Q ‘ > x 7 ) is obtained. 
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In case (a), the sequence st of states identifies n distinct 
vertices Ba =e are eos of I(P) and returns to the initial 
vertex 2 If n s 2 then the generated inequalities oot ext) < 0, 
for i = 0, 2, ..., 2n-2, define a bounded (2,2)-polyhedron P having 
n vertices. P is contained in P so that I(P) ¢ I(P). Also, since 
each of these inequalities is satisfied for all x € I(P), we have 
I(Pc¢ I(P). Furthermore, each vertex of P is an integer point, 
so that P is integral. By Lemma 1-1, P = P'. If n = 2 then P' 
is the line segment having a and as endpoints. Only two 
inequalities are generated. They define the line containing P'. 

In case (b), the sequence st of states identifies n+l distinct 
vertices me ee eee, Re of I(P). Beyond vertex <— we discover 
an infinite number of boundary points of I(P). It follows that 

-P is unbounded. The sequence S of states then identifies n'+1 


2 


7 -2n' 
distinct vertices = K 5 cee, X 2 of I(P). Beyond vertex 


or we discover an infinite number of boundary points of I(P). 
The generated inequalities, eG) <0 for i=0, 2, ..., 2n 
and ot ts (xex®) <0 fori =0, -2, ..., -2n', define an unbounded 
(2,2)-polyhedron P having n+ n' + 1 vertices. An argument 
similar to the one given for case (a) establishes that P = P'. 

In the following we show that the integralization process 
terminates after a finite number of states have been generated. 

Consider the sequence st of + successors of an initial state 


(9°, eye We divide st into subsequences $5 vag Si eee, where 
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Aa for k = 0, 1, 2, ..., consists of the consecutive block 
of states having the same supporting corner oX as first component. 
Each supporting corner gX is defined by a system of two inequalities 


kx < a, According to Lemma 5-10, cg’)! has at most A, = | det c* 


+ 


vertices. It follows that the subsequence S. contains at most 


A. States. 
: + + + 
Now consider the subsequences So 7 Sy 2» Spo eee Associated 
with each subsequence is a supporting corner gk. Associated with 
corner go is an objective inequality drawn from the system of m 
inequalities Ax < b defining P. This inequality is identified by 


the subscript i, Consider the sequence jg» jy> jo» eee « There 


are two cases to consider: 


(i) The sequence contains a jy such that 4, = j. for some i < k. 


i 
More simply, the sequence contains a repeat. We assume that i, is 

the first repeat in the sequence. In this case, it can be shown 

that =" is a vertex of 1(Q*). It follows that, for some state 

aK, oO in ac 8 = —: If n = 0 then we have the degenerate 

case in which I(P) consists of the single integer point ae Otherwise, 
n> 2 (the case n = 1 is not possible since s- is necessarily distinct 
from <°) and we have a state, as in case (a) above, at which st terminates. 
Thus 7 is the final subsequence of st, Since each member of the 


sequence jg» 4y> jg; eee, jy, has one of m possible values and J, is the 


first repeat, it follows that k < nm. 
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(ii) The sequence Jo» dy> igs ++ contains a j, which 
identifies an objective inequality drawn from Ax < b whose boundary 
line contains a half-line edge of P. In this case, some state 

k 2n, , +, : : 
(Q°, x ) in Ss. is a terminal state, as in case (b) above. 
Thus ee is the final subsequence of st, The sequence ip» jj> ses. jy, 
contains no repeats, otherwise case (i) would apply. It follows 
that kem. 
- + + g tT 
In either case, the number of subsequences So ; Sy zcees SL 
+ 

in S is at mostm+1. Thus st is a finite sequence of states. 
A similar argument can be made for a sequence S of = successors. 
It follows that the integralization process terminates after a finite 
number of states have been generated. 


We conclude this chapter with an example illustrating the 


integralization process for a (2,2)-polyhedron. 


Example 6-2 Let P be defined by 


-1 0 x) 0 
s 

0) -1 X5 0 

-18 8 21 

-9 12 56 

4 5 45 
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We reduce P to the following intersection of integral half-planes, 


- ) 

1 xy 0 

s 

0 -1 Xo ) 
ee 10| — 
-3 4 18 

4 5 45 

8 3 68 


An initial state (or. =) is found to be, 


We generate a sequence of + successors of (Q°, re The steps 
involved are displayed by giving a tabulation of the resulting 

Za and sy vectors, followed by the sequence of supporting corners 
and labeled W matrices. In the tabulation, we circle the component 
of zt for which be is smallest. Inequalities defining P' are 


enclosed in boxes as they are generated. 
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The sequence of + successors terminates with (@ 


0 
=X. 


Fig. 6-3 illustrates P and H(I(P)). 


Yr 
YO 
WYYYYYNS 
YYGYYGGOLE 
WUGUGULUIN 


Fig. 6-3 


16 : 
, x ) since 
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contained in the non-negative orthant). Extending this algebraic 
characterization to include all totally integral polyhedra would 
be an interesting problem for future study. 

We have answered the second of the two questions above for 
all (n,k)-polyhedra, where k < 2. We began in Chapter 4 by showing 
that an arbitrary (n,k)-polyhedron can be transformed into a 
(k,k)-polyhedron by a transformation which preserves integralization. 
This enables us to focus on the problem of integralizing (k,k)-polyhedra 
without loss of generality. We then described the solution to the 
simple case in which k = l. 

In Chapter 5 we developed our main integralization results for 
the case k = 2. We considered the problem of integralizing a (2,2)-corner 
polyhedron. Using a transformation which preserves the hull-forming 
operation, we showed that vectors normal to the surface of the sought 
after integral polyhedron are necessarily atoms of a certain partially 
ordered set. We introduced a two-dimensional generalization of 
the division theorem for integers, and showed how this is used 
to generate the finite set of all such atoms. We then described 
a second level division process involving these atoms which 
generates all vertices and normal vectors of the integral polyhdron, 
and thus its defining system of inequalities. 

In Chapter 6 we considered the integralization problem for 


(2,2)-polyhedra in general. We showed how a sequence of supporting 
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(2,2)-corner polyhedra is generated and how the integralization 
method of Chapter 5 is applied to each, thereby producing all 
vertices and normal vectors of the sought after integral 
polyhedron, and its defining system of inequalities. 

The integralization process developed in Chapters 5 and 6 is 
number theoretic in nature. It makes use of simple arithmetic 
operations without recourse to exhaustive enumerative methods 
or search procedures. 

It had been our goal to generalize this integralization 
method to make it apply to (k,k)-polyhedra, where k is arbitrary. 
However, we encountered some difficulty in attempting to generalize 
to k = 3 which we now describe. 

Let P be a (3,3)-corner polyhedron defined by Ax < b, where 
A is an integer 3 y 3 matrix with non-zero determinant and b is 
an integer 3-vector. We may integralize each of the three 
(3,1)-polyhedra (half-spaces) forming P. This gives us a 
(3,3)-corner polyhedron pi defined by a system aly < nie of 


three inequalities for which 
i 1:,1 
ged [Ag] = ged IA, | bg] 


for all face l-subsets S. We may then integralize each of the three 


(3,2)-polyhedra defined by two inequalities drawn from atx < 5". 


Thus gives us a (3,3)-polyhedron p* defined by the system APs < b 


of three or more inequalities for which 
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21 
ged [ac] = gcd [A, tbe] 


for all face 2=subsets S. Thus far we have used one and two-dimensional 
integralization methods to obtain ee. The remaining integralization 
step is three~dimensional in nature and this is where the difficulty 
arises. 

To illustrate the difficulty in its simplest form, suppose P 


is a (3,3)-corner polyhedron defined by Ax < b, where 


ecd [Ag] = ged [A, ‘bo ] 


S 


for all face l=-subsets and 2-subsets S. Then P = pt = P’, We regard 


A as a non-singular linear transformation, and let I(P)* be the image 
of I(P) under A. Normal vectors and vertices of L(P)* are defined 
using obvious generalizations of definitions of their two-dimensional 
counterparts, We let A be the cofactor matrix of A. As before | 
it can be shown that normal vectors of 1(p)* are necessarily members 
of M(A)*. However, it is unfortunately no longer true that normal 
vectors of 1(P)* are necessarily atoms of the partially ordered set 


(A )*, <). The following example serves as a counter-example. 


Example 7-1 Let P be defined by 


-4 i 1 x 3 
3 W5 4 Xo < 3 7 
3 -3 2 Ky 3 
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It is easy to verify that 


gcd [Ag] = gcd [A, tba] 


) 


for all l-subsets and 2-subsets S. The cofactor matrix A is given by 


2 6 6 
re 5 -ll -9/, 
9 19 17 


* 
Fig. 7-1 below illustrates H(I(P) ), determined by actual construction. 


(03,3) (2,33) 


13 
4 


a 


NORMAL 
Vector 
(4,2,2) 


(3, 3,-3) 
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The vector y = (4,2,2) is a normal vector of (p)* at vertices 

(3,1,;3), (3,3,1) and (2,3,3). It is easily verified that 

Y ce M(A)*. However it is also easily verified that (4,0,0) and 

(0,2,2) are members of M(A)*. Thus the normal vector y is not 

an atom of (M( A)t, <). 

O 
Since normal vectors are not necessarily atoms, our 

two-dimensional integralization method does not generalize 

directly. In order to generalize the method it would first 

be necessary to identify some finite subset of MCA )* which 

contains all possible normal vectors of (P)". Then a method 

for generating this finite subset of M(A)* would be required, 

along with a process for using this finite subset to generate 

all normal vectors and vertices of 1(P)*. In the two-dimensional 

case, it was the finite set of atoms of (M(A )*, <) that satisfied 

our first prerequisite. In the three-dimensional case, the 

appropriate finite subset of M( A) may possibly be those members 

of M( A)* which occupy the first two ranks of the partially ordered 

set (MC Ay <). This is merely a conjecture which seems to be 

a reasonable inductive inference to make from the two-dimensional 

case, once the set of atoms of (M(A yt, <) are found to be insufficient. 
Our two-dimensional integralization method gives us reason to 


believe that there exist number-theoretic methods for integralizing 
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(k,k)-polyhedra for arbitrary k. Our experience in seeking a 
generalization leads us to believe that inventing these methods 
is apt to be difficult. The problem remains a fascinating one 
for us, and one whose solution we believe would be of great value 


in the field of combinatorial optimization. 


3. 
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